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ON THK ORIGIN OF THE SOLAR SYSTEM 
By GERARD P. KUIPER 
YERKES OBSERVATORY, UNIVERSITY OF CHICAGO 


Read before the Academy, October 12, 1950 


Statement of the Problem.—-A satisfactory theory of the origin of the solar 
system must account for the presence and the properties of the planets and 
the smaller bodies surrounding the sun, and preferably, but not necessarily, 
for the dynamical properties of the sun also. This means that we shall be 
concerned with the following bodies and properties: 

(1) Nine Planets: The orbits around the sun are nearly coplanar and 
nearly circular; the largest inclinations and eccentricities are those of the 


outermost and innermost planets, Pluto and Mercury. The motions 
around the sun are all in the same sense (direct, by definition). The rota- 
tion of the planets is also direct, with obliquities less than 30°, except for 
Uranus where the angle has the exceptionally high value of 97°. The dis- 
tances of the planets from the sun exhibit some degree of regularity (Bode’s 
law). The masses of the four inner planets are roughly 10°-*© and the 
densities 4.1-5.5 cgs.; while the four Jovian planets have masses some 
hundred times larger and densities between 0.7 and 2.5 egs. 

(2) Thirty Known Satellites: The satellite systems vary from the 
beautifully regular case of Uranus to a completely irregular system like 
Neptune and the abnormal Earth-Moon system. Partly regular and 
partly irregular systems are those of Jupiter and Saturn. ‘“‘Regularity’’ is 
measured by low inclinations with respect to the planetary equator; small 
orbital eccentricities, e; direct motion with respect to the planetary rota- 
tion; and some degree of regularity in the mean distances, a, to the planet. 

(3) The Asteroids and Meteorites: The a-values of the 1168 largest 
asteroids! are roughly distributed according to a = 2.89 + 0.24 (p. e.) 
astr. units, though irregularities and fine structure occur in the distribution 
(Kirkwood Gaps) caused by the perturbations of Jupiter. Some 30,000 
asteroids are accessible to telescopic observation and the indications are 
that their a distribution is not very different from that of the largest mem- 
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bers. Meteorites strike the Earth after having arrived in orbits showing a 
distinct similarity to those of asteroids (i ~ 10°, e& 0.2). 

(4) The Comets, Both Periodic and Non-periodic, with Their Associated 
Meteor Showers: The i-values are nearly random except for the short- 
period comets and meteors which show a preference for small 7 and direct 
motion. The e-values are moderate for the short-period comets and almost 
unity for the remainder. Since the long-period comets with e < 0.9 cannot 
approach the Earth, we have no direct knowledge of them. 

The sun has an angular momentum of only about 0.3% of that possessed 
by the entire planetary system. The smallness of this ratio requires an 
explanation, but the process controlling the origin of the planetary system 
is not necessarily the sole cause of it. 

Analogy with Binary Stars.—The solar system is not the only multiple 
system known to astronomy. Simple calculation shows that planetary 
systems other than our own would be inaccessible to observation. But 
companions massive or luminous enough to be observable are found in great 
numbers. <A detailed study of the nearest stars shows that at least half of 
these are in fact double or multiple systems. 

The mean separations a in binaries are roughly given by? an error curve in 
log a having a dispersion in “log a of \/2: 


Filog a) Ce [log a~1.27)/2]? (1) 


The median value of a is therefore about 1S astr. units, a value remarkably 
close to the position of the most massive planets in the solar system. 

The mass ratios in binary systems of approximately the solar mass have a 
distribution roughly given by*® 


Moe 
if = constant (vatia for 0.3 < < 1.0) (2) 


my, my 


m 
The data for 0 < ~* < 0.3 are still too sporadic to give a trustworthy check 
my 


on equation (2) in that range. But several mass ratios down to about 0.02 
have been found and equation (2) appears to give the correct order of 
magnitude of the relative frequency of such cases. The absolute frequency 
of binaries with large mass ratios is found by putting the constant in equa- 
tion (2) at about 0.7, 1.e., a quantity of the order of 1. We shall, therefore, 
obtain a tentative estimate of the occurrence of small companions by using 


= 1, (2’) 
m, 
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Anticipating the coriclusion made below that the formation of the plane- 
tary system was a special case of the very general process of binary forma- 
tion, we may estimate the probability of this event by means of equation 
(2’). As the appropriate mass m. we may use either the mass of Jupiter, 
10°-*©, or the hypothetical mass of the entire solar nebula from which all 
the planets formed, 10°-'©. The use of the second quantity requires a 
further estimate of the fraction of cases in which a secondary nebula of mass 
10-'© was spread out and formed several small bodies instead of a single 
stellar companion. This fraction, obviously <1, may as to order of magni- 
tude be put at 10°-'. The probability that a main-sequence star like the 
sun is attended by a planetary system may therefore be estimated to be of 
the order of 10? or 10 *. On the basis of these developments the event 
would be enormously more probable than on the basis of the Tidal Theory 
or any other theory requiring a stellar collision or near collision; the latter 
would make the probability per star of the order of 10~'°'*, or 10* times 
smaller than according to the process envisaged here. 

The problem of the origin of binaries is considered in a separate article; 
in it a provisional interpretation is given of equations (1) and (2). It is 
shown that virtually all stars are likely to have originated as components of 
double and multiple systems. This follows from the result that the proba- 
bility of a proto-star having a small-enough angular momentum to allow 
the formation of a single, rotationally stable, star is less than 10~%. Many 
of the original double and multiple stars must have dissolved iater due to 
internal gravitational perturbations and effects of passing stars, so that at 
present the fraction of apparently single starsis substantial, roughly '/, to'/s. 

The analogy between binaries and solar system, suggested by the simi- 
larity of internal dimensions, will be complete if it can be shown that a 
small pre-stellar cloud moving about the primitive sun may form either a 
single stellar companion or a family of planets, depending on the mass dis- 
tribution within this secondary cloud. 

The Solar Nebula. The composition of the planets, particularly of the 
terrestrial planets, differs appreciably from that of the sun. It is now 
realized that the composition of the sun is typical for that of stars and even 
interstellar matter; it may therefore be called the cosmic composition. 
The cloud from which the planets formed will almost certainly have been of 
cosmic composition. A comparison of the planetary composition with the 


cosmic composition will therefore give a lower limit to the mass of the 
nebula from which the planets condensed and will contain many clues to 
the condensation process itself. In this manner the mass of the solar 
nebula may be estimated to be in excess of 0.01 © and probably more nearly 
0.1©. The realization that the pre-planetary cloud (or solar nebula) had 
this comparatively large mass has caused a revision in our concept of the 
possible processes by which planetary condensation may have taken place. 
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It appears now that gravitational actioa of the nebula on itself can have 
been important. It is seen below that precisely this property has made the 
formation of planets possible at all. 

We shall now sketch the principal processes leading up to the formation 
of the planets, their masses, orbits, rotations and satellites. A more com- 
plete account is given in two papers now in press.° 

Celestial mechanics indicates that the planets have not altered their dis- 
tances from the sun by large factors.® (Later we shall find that the altera- 
tions were, in fact, remarkably small.) The solar nebula must therefore 
have extended at least from Mercury to Pluto. It appears that the density 
distribution in this nebula need not be assumed, but can actually be de- 
rived from the present planetary distances; but as a first step we may sim- 
ply assume that the total mass of the nebula is 0.10, as is suggested by the 
present planetary masses and compositions. This assumption had already 
been made by v. Weizsaécker? who found that turbulent dissipation will 
cause the solar nebula to dissolve with a half-life of about 107 years. The 
dissipation is caused by the different angular velocities of the different 
parts of the nebula. The time available for the formation of condensation 
products is therefore limited. 

Another process, working in the opposite sense, is the slow collapse of the 
nebula parallel to the resultant axis of rotation (the 2-coordinate). The 
initial stages of the collapse are fast, everywhere comparable to the orbital 
period of rotation around the sun. The later stages are much slower, 10‘ 
10° years in duration, and governed by the rate at which excess heat, 
caused by contraction, can be radiated away. Therefore, a maximum of 
density will result, some 10*-10° years after the solar nebula got first in 
place; it is expected that the planets segregated close to or somewhat before 
that epoch. 

Condensation products will form wherever the vapor is supersaturated.* 
Unless some large-scale coalescing property within the solar nebula can be 
found we must anticipate the formation of enormous numbers of small 
condensations (a process somewhat comparable with the formation of rain- 
drops in a cloud), each growing slowly with time. The coalescing property 
needed to account for a small number of large planets was assumed by 
Weizsicker to be a regular turbulence pattern in the solar nebula. Weiz- 
sicker assumed a geometrical pattern of primary vortices, separated by 
rings of secondary “‘roller-bearing’’ vortices; he further assumed that con- 
densation would progress most rapidly along the roller-bearing circles. 
This was assumed to lead to masses of sufficient prominence to insure their 
stability against subsequent destruction and hence their subsequent 
growth by accretion of neighboring material. Ter Haar® '° elaborated 
Weizsacker’s picture and refined the discussion of the condensation process. 
Chandrasekhar and Ter Haar'® showed that the regular turbulence pattern 


. 
t 
q 
tease = 


Vor, 37, 1951 ASTRONOMY: G. P. KUIPER 


envisaged by Weizsécker would not arise but should be replaced by the 
Kolmogoroff spectrum of turbulence previously derived for non-rotating 
media. The writer thereupon showed® that secondary eddies were less 
suited to foster condensation than the primary eddies of the Kolmogoroff 
spectrum, because any “‘roller-bearing’’ eddies that may form when large 
clouds collide would be short-lived and of distinctly higher temperature, but 
only slightly higher density, than the primary eddies. He further showed 
that even the primary vortices had comparatively short lifetimes, of the 
order of 10? years rather than 10° years as is necessary for the formation of 
appreciable condensations. The repeated formation and dissolution of 
primary vortices in more or less random positions woud lead to innumerable 
small condensations throughout the solar nebula instead of a few large ones. 
This obstacle can be overcome only if it is assumed that gravitational action 
intervenes and keeps some of the large vortices together long enough for 
massive condensations to form. This requires that these vortices have a 
density large enough to protect them from tidal disruption by the sun, i.e., 
their density must exceed the critical value, pz, at which the internal attrac- 
tion within the cloud and the tidal force by the sun are equal. This critical 
density is called the Roche density; it is intimately related to the well- 
known Roche limit, used, e.g., in the discussion of the stability of satellites 
and the Ring of Saturn. The Roche density varies inversely as the cube of 
the distance from the sun and is given by 


3) 


= 68 IMo 
Pr = Po = = 


‘ne? 
mas GP? 


a 


in which 8 is a quantity of the order of unity, a the distance to the sun and P 
the corresponding period of revolution around the sun; Ro and ‘W/o are the 
radius and mass of the sun. The value of p, is about 10~° egs. near Mer- 
cury, 10°° near the Earth, near Jupiter and near Neptune. 
(The precise values depend not only on the mass and the distance to the sun 
but also on the molecular and turbulence velocities;> the range of uncer- 
tainty caused by the latter is covered by the approximate nature of 8.) 
Condensation of masses of planetary size could not have occurred unless 
the mean density of at least some vortices or clouds in the solar nebula ex- 
ceeded pg. On the other hand, the mean density of these clouds could not 
have exceeded p, by a large factor either. Such an excess in density would 
have caused the domain of the cloud to extend beyond its initial boundary 
and therefore the adjacent parts of the solar nebula would automatically 
have become part of the cloud until the mean density was lowered to a 
value of the order of p, itself. 

The conclusion to which we have come, that the segregation of pre- 
planetary masses took place essentially at the Roche density, appears to be 
capable of empirical verification. It can be shown® that the stability of 
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two small spherical masses in contact, each at the Roche density, is not 
appreciably affected by the presence of the remainder of the solar nebula 
except near its outer edge. Therefore if m and r are the masses and the 
radii of the small spheres, we find from (3): 


m 2rB\3 x 
= ) . (4) 
Mo d 
If we introduce A = 2r = ad, — ad, the distance between the centers of the 
masses whose mean distance from the sun is a, we have 


m (y. (4’) 
Mo a 

If A should in some cases be as large as a3 or a/2, as is indicated by the 
present planetary separations, we must allow for the fact that the flat.disk- 
shaped solar nebula would not have filled such large spheres except for an 
equatorial slice occupying the fraction f of the order of perhaps 107'. 
Further, a comparison of the cosmic abundances with the present chemical 
composition of the planets shows that no more than a fraction g of the 
original pre-planetary clouds has been retained, g being about 107! for 
Jupiter and 10~* for the Earth. The geometric mean of g for all planets 
may be roughly 10~*. Therefore, if 7» represents the present planetary 


masses, we should expect 

A\3 

Mr _ | 
Mo a 


while empirically we find approximately® '! 


A\3 
= 10 (5) 
Mo a 


It may be shown® that fis between 10 'and 10~*. The remarkable agree- 
ment between the predicted relation (4”) and the empirical relation (5) 
shows that the segregation of planetary masses did take place essentially at the 
Roche density. The closeness of the agreement is appreciated by remem- 
bering that the Roche density varies by the factor 500,000 between Mercury 
and Neptune and the \/, values are proportional to this quantity. 

It appears that the comparison between equations (4”) and (5) may even 
be refined to take account of the variation with J/, of gin (4"). For the 
pair Jupiter-Saturn the value of f is about 107 ''*, while g is roughly 107! 
If 8 = 1 we should for this pair expect the coefficient in (5) to be 1077", 
This is precisely what is found.'! For the pair Earth-Venus, f = 107}, 
g = 10-3, roughly, and hence the expected coefficient in (5), 10-4. This, 
again, is confirmed." It is significant that the empirical relationship (5) 
holds for both planets and satellites; we return to this comparison later. 
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These examples show that the systematic deviation'! from equation (5) 
is caused by the systematic change of g with Mp. A second result," that 
pairs with very unequal masses are systematically somewhat displaced 
with respect to pairs of nearly equal masses, has not yet been interpreted. 

With the density in the plane of symmetry determined, the density dis- 
tribution at right angles to this plane (z-direction) may be computed. 
Weizsacker and Ter Haar in their discussions of the hydrodynamics of the 
solar nebula assumed that the z-component of the force of gravity is given 
by the z-component of the solar attraction. With the equatorial density as 
high as pz the z-force near the plane of symmetry will in fact be predomi- 
nantly caused by the attraction of the nebula itself, in much the same way 
as the z-force in the galaxy near the galactic plane is determined by the 


i 
-6 -2 Log Mp, 
FIGURE 1 


Relation of computed mass of proto-planet, Mpp, to pres- 
ent observed mass of planet, Mp. Unit of mass is Sun. 


density distribution of the surrounding parts of the galaxy and not by the 
galactic nucleus. The ratio of the two z-forces, nebula/sun, is 248%, or 
about 25. The solar nebula is flatter than considered by previous investi- 
gators,” * '° by this same factor. The resulting large density increase in 
the new model allows the operation of gravitational instability. It is this 
mechanism which amplifies initial density fluctuations caused by turbu- 
lence’ and is responsible for the formation of large clouds (proto-planets) 
which are gravitationally stable. These stable clouds exist long enough 
for the condensation processes to run their course and for large masses to 
form instead of billions of small, unrelated drops or flakes. 

The density distribution in the z-coordinate may be computed from an 
equation ‘resembling the barometric formula, the field of force being known 


z 

° 
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and also, approximately, the molecular velocities. The latter follow from 
the distribution in temperature with distance from the sun, which in turn 
may be estimated when allowance is made for the opacity of the nebula. 
The condition that the proto-planets segregated at the Roche density, 
together with the density gradients in z just found, fixes completely the 
mass distribution for each proto-planet, both in the xy plane and the z- 
direction. We need only to add the dimensions in the xy plane for a com- 
plete description of the proto-planets, including their masses. The diam- 
eters of the proto-planets, at their maximum extent, are found by assuming 
that the proto-planets grew until they touched their neighbors on each 
side and by the further assumption that the planets formed at their present 


TABLE 1 
LoG a a MAX. ~ M(PP) 
PLANET a BOUNDARY ALOG 4 BOUNDARY 2R(PP) M(PP) M(P) 
Mercury —( 330 
(aph.) —0.24 0.58 

Venus —(). 141 0:47 0.27 0.0037 3.18 
—0.07 0.85 

Farth 0000 0.16 0.388 0.0088 3.10 
+009 1.23 

Mars +0. 183 0.18: 0.68: 0.0045: 4.14: 
+0). 27: 1 86: 

Asteroids +O 45: 0). 28: 
+0). 55: 3.55: 

Jupiter +0. 716 0.30: 3.5 0.012 1. 0 
+0. 85 

Saturn £0,979 0.80 7 0. 014 1.69 
41.15 14 

Uranus +1283 0.23 10 0 0072 2.22 
+1.38 24 

Neptune +1.478 0.20: 14 0.0057: 2.04 
+1.58: 38 

Pluto +1. 597 


NorkE: Unitofaisa.u. PP = proto-planet; ? = planet. Cf. figure 1. 


distances to the sun. The resulting partitioning of space is shown in table 
1; it was done in a logarithmic plot because of the approximate validity of 
Bode’s law and the fact that all relevant functions have an exponential run 
with a. The hypothesis that the proto-planets at their maximum size 
touched each other, accounts for the fact that at that point further growth 
was arrested; nearly all of the solar nebula would thus have been ‘swept 
clean.”’ At the same time a mechanism is provided for the explanation of 
retrograde satellites in the outer regions of the “sphere of action’’ of such 
planets as did not lose thereafter too great a fraction of their original mass 


(Jupiter and Saturn). 
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The preceding arguments fix the masses of the proto-planets, J/ pp, 
uniquely. They are found in table | and plotted against the present plane- 
tary masses, \/,, in figure 1. The closeness of the relation is very remark- 
able. The largest deviation is that of Mars, amounting to 10°? or a factor 
1.6in Mpp. This may be interpreted to mean that the diameter of proto- 
Mars was in fact only 0.S of that assumed; this difference is not serious 
since the outer boundary of the Mars zone is ill-defined. 

The curve relating \/,, to \/, has been tentatively extended beyond the 
region of the planets, as follows. For any value of 17, we must have MW, < 
Mpp. The limiting line so defined is dotted in figure 1; very probably this 
line is an asymptote. Further, it has been shown® that there is a theo- 
retical lower limit to >», given by log M,p = —4.8. All points must 
necessarily lie above this limit. It is assumed that Ceres and a few other 
asteroids are true planetary condensations. Their 1/, values are about 
10-"©. The shaded region in figure | shows the most probable location of 
these objects. 

The closeness of the relation in figure | shows two things: (a) The position 
of Uranus and Neptune is intermediate between that of the Earth and 
Venus on the one hand and that of Jupiter and Saturn on the other. ‘This 
shows that the relation in figure | is not strongly affected by the distance 
to the sun, though an extra factor of two or three in the present masses of 
the Jovian planets would hardly be discovered in this manner. This 
means, e.g., that the Earth is different from Jupiter not because it originated 
closer to the sun but primarily because proto-Earth was less massive than 
proto-Jupiter. The existence of a planet like Pluto outside Neptune be- 
comes at once “‘acceptable.”” (b) The distances of the planets to the sun 
cannot have greatly changed during the last 3.10° years, at least not in an 
irregular fashion. If they had, the relation in figure | could not possibly 
have been close since |/,p depends in a sensitive way on the present plane- 
tary distances. This conclusion goes much beyond the powers of celestial 
mechanics which for such extremely long intervals has made only the most 
general predictions.*® 

The masses of the proto-planets are about as large (within 20 or 30%) as 
might have been expected from the present planetary compositions. This 
conclusion is an important check on the theory presented. It is the first 
time that planetary masses have been predicted with any degree of accu- 
racy. 

By allowing 0.003 M for each of the proto-planets which formed Mer- 
cury, the asteroids and Pluto we find from table | the total mass of that part 
of the solar nebula that was responsible for the formation of the planets. 
It is found to be 


> = 0.060 Mo. (6) 
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The total mass of the solar nebula must have been somewhat (but probably 
not appreciably) larger and may be put in round figures at 107! Mo. 

The Bode Law. Planets and Binaries.—We must now consider the fac- 
tors that determine the relative spacings, A/a, of the planets. Why did 
nine planets form and not 10? smaller ones? 

The mode of formation of the proto-planets gives a clue. If the model of 
one proto-planet per ring, valid for the large planets, is extended to smaller 
and smaller bodies, we note that a smaller and smaller fraction of the ring is 
occupied to the Roche density. It follows that the spacing of the planets 
at the distance a from the sun is determined by the ratio of the density in 
the solar nebula to the local Roche density. 


0.10 


loga +2 


FIGURE 2 


Mass distribution in solar nebula prior to formation of 
proto-planets (tentative ). 


In this manner the density distribution is found which existed in the solar 
nebula prior to the break-up into proto-planets. In order to express the 
result we take as reference a nebula having everywhere the Roche density 
in the xy plane. Such a nebula will have a given “‘surface density” of 
matter when projected on the xy plane. In fact, it is one in which the 
amount of matter per unit interval of log a is constant. In terms of this 
standard surface distribution, the pre-planetary nebula had the density 
0.05 from about '/2 to 2 astr. units, rising to 0.09 at 10 astr. units, and falling 
off to below 0.05 at 30 astr. units. It 1s illustrated in figure 2. 

From the theory of gravitational instability it follows that a lower limit 
exists below which stable proto-planets cannot form. This limit in the 
units of figure 2 is found to be 
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Co = 


in which both a and 8 are quantities of the order of unity. It appears there- 
fore that the solar nebula had only just enough surface density to cause 
formation of planets. A reduction by a factor of two or three (or more) 
would have been fatal to this process; billions of comet-like bodies would 
have formed instead. This result may also be interpreted to mean that 
the hydrogen content of the solar nebula could not have been much less 
than supposed; a reduction by a factor three over the adopted cosmic 
composition would have caused a nearly equal reduction in mass, and 
would again have prevented planets to form. It is possible that the 
hydrogen content was actually somewhat larger than assumed. 

If, on the other hand, the surface density of matter over the xy plane had 
been larger than found in figure 2 by a factor of three or more, no ordinary 
planets would have formed either but instead one or more stellar compan- 
ions. This results from the circumstance that a rather small increase in C 
(see Fig. 2) causes a large increase in the mass of the proto-planet, which 
in turn causes a still greater increase in the mass of the final body (Fig. 1). 

Three conclusions follow: (a) the so-called Bode law is no real law but a 
consequence of the density distribution in a which the solar nebula hap- 
pened to have; the larger the density, in terms of the local Roche density, 
the larger the ratio between the radii of two consecutive orbits and the 
larger the masses of the proto-planets (and, finally, the larger the masses of 
the planets); (b) the process of formation of a planetary system is a special 
case of the almost universal process of binary-star formation; the crucial 
parameter is the surface density of the secondary nebula in terms of the 
local Roche density; (c) inspection of figure | suggests that the subsequent 
evolution of the proto-planets is not greatly influenced by the distance to 
the sun, but instead largely determined by the mass of the proto-planet it- 
self (once formed, the object is “‘on its own”). Space does not permit to 
discuss the theoretical interpretation of this conclusion; it appears to be 
connected with the rate at which the lighter gases can escape.° 

Conclusions.-We may now turn to the problems listed on page 1 and 
list the solutions now at hand or indicated. The common direction of 
revolution and the low relative orbital inclinations are accounted for by the 
flatness of the solar nebula. The internal viscosity of the nebula accounts 
for the near-circular orbits. The fact that both Mercury and Pluto are 
exceptions in their inclinations and eccentricities may be attributed to the 
absence of constraining action on the proto-planets formed on the fringes 
of the solar nebula. The direct rotation of the planets is attributed® to 
solar tidal friction on the proto-planets. The solar tidal force nearly 
equals the self-attraction for each of the proto-planets at their maximum 
extension; i.e., in the proper units Neptune is no farther away from the sun 


q 11 
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than Mercury and the tidal effects are equally large in both cases. Regard- 
less of any initial rotational motion a direct rotation will be forced upon 
the proto-planet, with a period equal to the orbital period. As is readily 
verified, this leads to an amount of angular momentum per unit mass some 
104 times greater than found on the present planets. Part of this is lost 
during the evaporation process of the proto-planets (the ejected molecules 
carry off more than the average amount of angular momentum per unit 
mass); while part of it is lost by continued solar tidal friction during the 
contraction process, The latter cause has a secular effect on the obliquities; 
it has been shown’ that they will increase some three or fivefold, from 
initial obliquities of the order of 3° (expected from the turbulent solar 
nebula, and consistent with the relative orbital inclinations) to the present 
values. The largest obliquity to which this process can lead is 90°; retro- 
grade rotation cannot arise by the processes considered. It is not clear 
why Uranus has passed the upper limit by 7°; possibly some extraneous 
object has moved through the solar system. The present periods of rota- 
tion have not yet been accounted for quantitatively. This appears to be a 
very complex problem, with physics, chemistry and dynamics all playing a 
role. We have here perhaps the most important potential source of in- 
formation still unused in the reconstruction of the planetary condensation 
processes, 

The regular satellites may be explained in a manner analogous to that 
found for the planets themselves.° Little progress has been made so far 
with their condensation processes, which should prove very instructive in 
view of the large density differences known to exist among the satellites. 
The retrograde satellites of Jupiter and Saturn have been interpreted® as 
having been caused by glancing collisions between the corresponding 
proto-planets. They were assumed to have been retained by these large 
planets only because these planets lost a much smaller fraction of their 
initial mass. It is possible, however, that capture has played a role instead. 
This requires further investigation. The asteroids were not formed in a 
region of low density in the solar nebula. In such a region no planets of 
any kind could have formed. Rather we must assume that the density was 
well above Cj of equation (7), but that the formation of a normal-size proto- 
planet was prevented by proto-Jupiter (mass = 0.0120). It can be 
shown that in the presence of strong perturbations a small proto-planet, of 
a given density close to the local Roche density, is more stable than a large 
one of the same density. The total number of small proto-planets esti- 
mated to have formed in the region between Mars and Jupiter is between 5 
and 10. They formed small planets, like Ceres (ef. figure 1 and accom- 
panying discussion). It is assumed that two of these collided sometime 
during the last 3.10" years, an event having a sufficiently large probability. 
Thereafter secondary collisions became increasingly frequent. The most 
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recent of these collisions account for the Hirayama families. In this 
manner thousands of asteroids were formed, being the largest of the frag- 
ments, as well as billions of meteorites. !* 

The outermost region of the solar nebula, from 3S to 50 astr. units (i.e., 
just outside proto-Neptune), must have had a surface density below the 
limit set by equation (7). The temperature must have been about 5-10°K. 
when the solar nebula was still in existence (before the proto-planets were 
full grown), and about 40°K. thereafter. Condensation products (ices of 
HO, NH3;, CH,, etc.) must have formed, and the flakes must have slowly 
collected and formed larger aggregates, estimated to range up to | km. or 
more in size. The total condensable mass is about 10** g., but not all of this 
could be collected. These condensations appear to account for the comets, 
in size,'* number'* and composition.'* 

The planet Pluto, which sweeps through the whole zone from 30 to 50 
astr. units, is held responsible for having started the scattering of the 
comets throughout the solar system. Pluto's perturbations will have 
caused initial, near-circular, cometary orbits to become moderately ellipti- 
cal; thereupon stronger perturbations by Neptune and the other major 
planets will have scattered them even more broadly. As Oort'® and others 
have shown, the quantity which is spread nearly uniformly in both direc- 
tions is the quantity a~', the reciprocal of the semimajor axis (which is re- 
lated to the energy of the object). A certain fraction of the comets will be 
scattered in the region of very small a~! values, i.e., in the outer regions of 
the ‘‘sphere of action” of the sun. As Oort!* has shown, stellar perturba- 
tions will redistribute the orbital elements there, and in particular make the 
motion around the sun one of random orientation. Oort!® shows that the 
dynamical half-life of a comet in this outer region is about 10" years. The 
comets which we observe today were sent back to the inner regions of the 
solar system by small random stellar perturbations. The above views are 
an adaptation of Oort’s'? dynamical analysis; but we differ in our hypo- 
thesis as to the region where the comets originated. Oort'® assumes that 
they were formed between Mars and Jupiter, in association with the origin 
of asteroids. The composition of the comets indicates condensation at a 
very much lower temperature, around 10°K., consistent with the region of 
origin proposed here. The evaporation and subsequent complete disinte- 
gration of comets into the minute particles which cause meteors and the 
Zodiacal Light is also understandable from their formation outside Neptune. 
Asteroidal bodies would be expected to remain intact or possibly break up 
into a few large fragments. 

The theory described here does not depend on any specific ad hoc assump- 
tions. Certain assumptions which were made at the outset, e.g., that the 
planetary distances have not changed appreciably or that the solar nebula 
was approximately of cosmic composition, appeared capable of verification 
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afterwards. One assumption, that the sun was already formed as a star 
and of a luminosity approximately equal to that found today, requires 
further study.'® Certain investigations on the contraction and condensation 
process of the proto-planets need still be made, including the analysis of 
solar tidal friction on these composite structures. Finally, the cause of the 
small solar rotation must be cleared up; it is undoubtedly connected with 
the larger problem of why nearly all G-type dwarf stars, in single and in 
binary systems, have such slow rotations. It is felt, therefore, that this 
problem is not necessarily a part of a theory on the origin of the solar sys- 
tem. 

The probability of a star being attended by a planetary system was esti- 
mated to be between 10~* and 10~*. The total mass of the galaxy is about 
2.10''O; while the average stellar mass is about 0.5©. From these figures 
the total number of planetary systems in the galaxy is estimated to be of 
the order of 10". One can only speculate on the possible forms of life 
which may have developed on these many unknown worlds. 


' Those with g < 10.9 in the 1947 Catalogue of Minor Planets. The quantity g is 
the magnitude of the asteroid reduced to unit distance from the Earth and unit distance 
from the Sun; it depends chiefly on the size of the body since the albedos are probably 
not very different. 
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PROXIMA CENTAURI AS A FLARE STAR 
By HARLOW SHAPLEY 
HARVARD COLLEGE OBSERVATORY 
Communicated November 11, 1950 


1. The dwarf red companion of the bright double star Alpha Centauri 
was detected by Innes in the course of an examination of plates made 
with the Franklin Adams camera at the Union Observatory, Johannesburg, 
in 1915.1. Although the object was two degrees distant from Alpha 
Centauri, its similar motion, and subsequently its comparable parallax, 
clearly showed the physical association of the three objects. The some- 
what larger proper motion and the slightly larger parallax, derived by some 
of the observers, have led to the conclusion that Proxima Centauri is 
actually a little nearer than Alpha Centauri, and is therefore the nearest 
of all stars so tar discovered; hence the appropriateness of its name. The 
difference in the parallax, however, scarcely exceeds the probable error of 
the determinations, and van de Kamp gives in his listing of the stars nearer 
than five parsecs* identical values for the parallax, 07761, and for the 
annual motion of Alpha, 3768, and of Proxima, 3785. 

2. It became clear some years ago that such low-luminosity red dwarfs, 
at the faintest end of the main sequence, are not simple and quietly perish- 
ing stars. Castor C, for example, was found at Mount Wilson to be a 
close pair of dwarf Class M stars. It is an eclipsing variable with a period 
of O.S14 days. Later, in the course of his parallax and proper motion work, 
van Maanen found that Lalande 2125SB and Ross S82 had each shown one 
abnormal brightening among the many images appearing on the Mount 
Wilson plates.* 

In 1949, Luyten’s remarkable double, L 726-8, discovered by him as a 
large proper motion star on the Harvard Bruce plates, was found to be a 
light variable. Not only is it one of the nearest of all stars, and among the 
dimmest stellar radiators yet on record, but also the fainter component 
of the pair showed brief flares on two occasions. The star has been fully 
discussed by Luyten.! In their concurrent work on the spectra, Joy and 
Humason? also found a flare, especially strong in blue light. Subsequently 
569 Harvard plates yielded seven additional flares,* one of them corre- 
sponding to a momentary increase of forty times in the light intensity. The 
results for L 726-S not only gave support to van Maanen’s earlier observa- 
tions, but attracted general attention to the dwarf red stars. There are not 
many of them on record, and they cannot be seen very deep in space; 
but here were three affected by temporary flashes-~a novelty in the photom- 
etry of stars. The spectrum of Luyten’s star showed unusual charac- 
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teristics, including bright calcium and hydrogen lines.° Obviously these 
red stars are not cooling off monotonously. 

The list was increased to four by the interesting photoelectric discovery 
and observation by Gordon and Kron of an hour-long flaring-up of the 
M dwarf, BD + 20°2465, which confirmed earlier evidence for the other 
flare stars that the duration of the spurt and recovery in brightness is in- 
deed of short duration.’ 

3. Before any of the foregoing information on flare stars had become 
known, nearly a hundred Harvard photographs of the field of Proxima 
Centauri had been examined in search of variability of the continuous 
sort that affects the dwarf red variable of the Castor system. The results 
were negative. The constant photographic magnitude was established, 
with the aid of the most convenient magnitude sequence, as 13.4, with 
small dispersion around that value except for one very discordant measure 
which was ignored. When Thackeray recently published a discussion 
of the variable spectrum of Proxima* and reported also that the observers 
at the Cape Observatory had examined 58 plates of the star, taken for 
parallax during the past two decades, and found no certain evidence of 
flares, the earlier Harvard measures were reexamined and the one most 
discordant observation was recognized as a probable flare. Four measures 
of 13.2, 18.2, 13.0, and 13.2 possibly also indicate flare phenomena, and are 
included in the later tabulation. The plates had been made with the 
10-inch Metcalf on a neighboring Milky Way Field, No. 171, and Proxima 
is in a rather unsatisfactory position on the photographs. 

An examination of several series in the Harvard plate collection has now 
been undertaken, ‘The MF plates were searched by Miss C. D. Boyd and 
the other series by Mrs. V. McK. Nail. Without going further back 
than 1925, five hundred additional plates were found on which the image 
of Proxima Centauri appears. The following tabulation gives the plate 
series, name of telescope, the average exposure time, the number of plates 
examined, and the number of times that Proxima Centauri was found to be 
of magnitude 12.8 or brighter. 


MF Metcalf triplet 45 min, 93 plates 1 flare 
B Bache doublet 45 80 6 
RB Ross-Fecker lens a0 170 15 
AM Cooke lens. 90 211 19 
AX Tessar lens 120 35 7 
A Bruce doublet 120 3 0 


Total 592 48 


It appears that on eight per cent of the plates (six per cent of the exposure 
time) the star was brighter than average by a distinctly measurable amount, 
making it the most active flare star on record. By going to the older 
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Harvard plates the list of recorded flares could be substantially increased. 
Possibly such a study should be made if it later appears that there is an 
indication of periodicity in the outbursts, but at present it is more important 
to follow some of these flares carefully with photoelectric photometers or 
short exposure photographs in order to learn more about the nature of the 
outbursts than can be obtained from these long-exposure Harvard plates. 
Proxima has now been put on the program of the photoelectric photometer 
at the Harvard station in South Africa. 

For all series the exposures are so long that brief sharp maxima are 
much smoothed out. The brightest recorded maxima are 12.4, indicating 
amplitudes of at least one magnitude. The AX camera has a steep color 
curve and is not very reliable for work on such red stars. 

Oni three nights consecutive exposures were made throughout 5.5, 2.25, 
and 7 hours with the 10-inch Metcalf telescope. No measurable change of 
brightness was shown. In addition, on seven nights plates with different 
instruments are available. Four pairs occur at times of flares; they give 
magnitudes consistent with the presumption of outbursts of an hour or so. 

The Julian Day and decimal of a day for the flares on MF, B, RB, AM, 
and AX plates, with magnitude estimates by Mrs. Nail, are as follows: 

Jb AND GMT SERIES 3. JD AND GMT SERIES MAG JD AND GMT SERIES MAG 
2423996 534 AX 2.6 2428040.215 B 2.6 2431184.427 AM 12.7 

24244 827 AX 2.5 28276.553 RB y 31197.562 12. 

24286 694 AX 2.¢ 28283.542 B 2 31311.343 AM 12.6 

24294 . 693 AX 28289.417 AM 31330.264 AM 12 

25446. 285 AX y 28896.552 AM 2.; 31494.591 AM 12. 

26091 .372 MF 2.6 28994.495 AM 2.6 31584.310 B 12. 

26091 .422 RB 29051.408 AM 2 31614.414 AM 12 

26131 .331 MF 3. 29131.238 RB 2. 31620.287 AM 12.8 

26181 MF 3.2 29834.345 AM 2.4: 31873.323 AM  12.- 

26189. 286 RB ‘ 29850 .233 RB 2. 31887.529 AX 12 

26215. 220 MF 3. 30134.463 AM 2. 31921.424 AM 12 

26812 .531 RB 30474.446 RB 31935.410 RB 12. 

27519. 505 AX 2.4 30491.354 AM 2.6 32033.290 AM = 12 

27525 635 MF 3. 30497 .277 2. 32324.277 12.6 

27630 . 224 RB 2.6 30515.448 RB 2.6 32387 . 291 AM $12.6 

27958. 469 RB 30578. 280 AM ‘ 32406.227 AM = 12 

28010. 320 RB 2.8 80850. 451 RB y 33033.577 RB 12, 

30900.320 RB 


4. A sixth star of spectral Class M, which has the spectrum charac- 
teristics of SS Cygni, may be assignable to this small group of flare stars. 
It is AE Aquarii, which has been studied at various observatories and is 
now on the program of the American Association of Variable Star Ob- 
servers.? Recent unpublished observations show rapid changes in the 
course of an hour or so. Like Proxima Centauri, it merits study with 


photoelectric photometers aud high speed spectrographs, 
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5. Dwarf red flare stars may become of considerable importance in 
considerations of stellar evolution. Something near one-half of the fifty 
nearest stars now known are dwarfs of Class M, with absolute photographic 
magnitudes fainter than +12. Of these, about twenty per cent are now 
known to be affected with sudden bursts of radiation, presumably emanat- 
ing, as Gordon and Kron, Whipple, Greenstein and others have pointed 
out, from relatively small areas which for a time explode with highly 
intense (blue) radiation, and quickly cool to the average conditions. 

The great contrast in surface brightness between the intense blue out- 
burst and the relatively cool surrounding surface greatly increases the total 
light of the red dwarf stars, whereas similar phenomena on the surface of 
the sun, or on stars of earlier spectral class, would not appreciably affect 
the total radiation. The total energy of an outburst on Proxima is equiva- 
lent to that of some of the solar flares currently observed; and although 
the solar flares do not measurably increase the magnitude of the sun, they 
do succeed in disturbing the earth’s ionosphere. 

6. Is it likely that measurable eruptive flaring disturbs the radiation of 
all the red dwarf stars: Possibly the phenomenon is universal, especially 
for red dwarfs with emission spectra, but the examination of 225 Harvard 
plates of the RB, RH and BM series on Barnard’s Star (spectrum M 6, 
absolute magnitude +13.4, vis.) reveals no variability. 

' Innes, R. T. A., Union Obs. Circ., 30 (1915), and 40 (1917). 

* van de Kamp, P., Pop. Ast., 48, 297 (1940). 

3 van Maanen, A., Mt. Wilson Contr., 630 (1940); Ap. J., 91, 503 (1940); Pub. Astr. 
Soc. Pac., 57, 216 (1945). 

‘Luyten, W. J., Ap. J., 109, 582 (1949). 

§ Joy, A. H., and Humason, M. L., Pub. Astr. Soc. Pac., 61, 133 (1949). 

® Harv. Ann. Card., 1010 (July 1949). 

7 Gordon, K. C., and Kron, G. E., Pub. Astr. Soc. Pac., 61, 210 (1949). 

§ Thackeray, A. D., Mon. Not., 110, 45 (1950). 

® Campbell, L., Pop. Ast., 55, 557 (1947). 

10 Greenstein, J., Pub. Astr. Soc. Pac., 62, 160 (1950). 
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THE THEORY OF MICRO-METEORITES.* PART II. IN 
HETEROTHERMAL ATMOSPHERES 


By Frep L. WHIPPLE 


HARVARD COLLEGE OBSERVATORY 


Communicated by Harlow Shapley, October 16, 1950 


Introduction._-In a previous paper! the writer presented a theory for the 
passage of a small meteoritic body through the earth's atmosphere without 
damage by melting or vaporization. The theory, which derives the maxi- 
mum ratio of mass to effective surface area for a micro-meteorite of a given 
initial velocity, is valid for an atmosphere of constant temperature and 
mean molecular weight but depends upon certain assumptions. The pres- 
ent paper extends the theory to more general atmospheric circumstances 
and also explores the tenability of the basic assumptions. For conven- 
ience, the previous notation is utilized while references and equations of the 
first paper are identified by asterisks. 

1. The Solution in an Atmosphere of Constant Temperature Gradient.— 
At altitudes in the neighborhood of 100 kilometers, most investigators at 
present tend to adopt an atmospheric temperature increase with altitude.’ 
Little is certainly known about the mean molecular weight at these alti- 
tudes but a constant value will probably be a fair assumption. The disso- 
ciation of O, probably occurs near 100 kilometers altitude but its effect on 
the density may be absorbed in the assumed temperature gradient for con- 
stant molecular weight, as may the slight effect arising from the decrease 
of gravity with height. 

It is, therefore, of great interest to repeat the solution for the maximum 
ratio of mass to surface area, applying to micro-meteorites in an atmos- 
phere with a constant temperature gradient. Let the atmospheric tem- 
perature, 7\,, vary with height, h, according to the relation 


Ta = 7, + vh, (1) 


where 7) is the temperature at a height, 4 = 0 and where v remains con- 


stant. 
The atmospheric density then can be expressed by the relation* 


h ~ (1 + q/v) 
p = pi (1 + (2) 


where p; is the density at h = 0. 
The quantity q is defined by 


ug/k, 
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where w is the mean molecular weight of the air, g is the attraction of earth 
gravity corrected for altitude and & is the Boltzmann Constant. The 
logarithmic density gradient, 4, used previously is equal to (q¢ + v)/T,. 

With our previously adopted assumptions that the drag coefficient, D, 
accommodation coeflicient, a, and emissivity, 8, are constants we may sub- 
stitute p from equation (2) for its isothermal value in equation (14*) and 
integrate the velocity differential to obtain 


ADT \ 
log (V/Ve) = — (: 4 *) (4) 


2mq cos Z 7) 


Equation (10*), giving the surface temperature of the micro-meteorite, 


aA p, V4 (1 + a/v) 


now becomes 


where we agai neglect ihe heat capacity. 

Upon eliminating the terms in /# from equations (4) and (5), we express 
the temperature of the micro meteorite in terms of the velocity, inde- 
pendently of height or density, as follows: 


V3 9 v/q os Z +q)/@ 
Ti ( ) | - log (V/ v.) | 


BBo Ap, D1, 


By equating to zero the derivative with respect to V of the right mem- 
ber of equation (6) we find that the maximum temperature of the meteor- 
oid, 7, occurs at the critical velocity, V., given by 

3q 


log (V./Va.) = — (7) 

This value of the velocity is now substituted into equation (6), A is 
equated to A, and expressed in terms of B by equation (2*) and the ratio of 
inass to surface area for the micro-meteorite is obtained in the form 


m + — Ty) +0 
(v+q)cosZ\8 


The reader will note that when vy = 0, eorresponding to an atmosphere 
of uniform temperature and molecular weight, equation (8) reduces to 
equation (18a*), since b = g/7\. Again the radius of a spherical meteor- 
oid of density, p, is s = (3m)/(p,B). 

For a temperature gradient of +4° per kilometer and V,, = 23.1 km./ 
sec., equation (8) yields a value of m/B some 10% greater than equation 
(18a*) for an identical atmospheric density and density gradient at the 
critical velocity in the isothermal case. Since these solutions are in such 
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close numerical agreement, we may trust equation (8) to give a rather pre- 
cise result for most problems of micro-meteorites. 

2. The Solution in the General Case.—-When the physical data justify 
a more precise theory for micro-meteorites than that given in the previous 
section, a numerical integration will be required generally. The atmos- 
pheric temperature will certainly vary with height in some complex fashion 
while its mean molecular weight may also vary appreciably. The drag 
coefficient and the emissivity may depend upon the velocity and tempera- 
ture of the micro-meteorite and possibly upon the air density. In addition 
the heat capacity, C,, may be appreciable and will vary with 7,. 

Before going on to the most general situation, let us investigate the error 
introduced by neglecting the heat capacity in the previous solutions. 
Beginning with equation (7*), the temperature differential can be related to 
the velocity differential by means of equation (13*), (14*) and (15*) so 
that the time, height and atmospheric density are explicitly eliminated. 
The temperature derivative with respect to the velocity is 

dT, w(T,4 — To! 
—*= —uV — 0) (9) 
dv V? log (V/V...) 
where the constants, # and w, are given by 
u = a/(CD) (10) 
and 
w = BBa/(bmC, cos Z). (11) 

Our previous assumption in neglecting the heat capacity was derived by 
equating the right member of equation (9) to zero. ‘The resulting equation 
(16*) 1s rewritten: 

T,4 — = —uV*(log V/V..)/w. (12) 

We may now derive a second approximation to the temperature by sub- 
stituting the derivative of equation (12) into equation (9). This operation 
leads to the expression 


To = —~ V* log (V/V) [ 


If we now differentiate equation (13) and set d7,/dV = O the result will 
be a second approximation to the velocity-temperature relation at maxi- 
mum temperature, 7). The equation is written more simply with the 
first order auxiliaries x and y where 


x = V,/(4wT,,?), (14) 


and 
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y = 1+ 3 log (V./Ve). (15) 


After dropping non-zero coefficients, we find the relation at maximum 
temperature to be 


y(1 — xy) = x(y — 1) (1 + y/3). (16) 
If we drop second-order terms in x and y, equation (16) simplifies to 


Hence, under the same conditions, 
log (V./V.) = (y — 1)/3, (18) 
and 
Vi = + y)/e. (19) 


Thus V,* is reduced by a first-order term (—y, which is positive). A 
substitution of the corrected value of V, into equation (13) to evaluate 
T.», however, leaves no first-order correction terms. The original value 
of 7,4 is unchanged except for second and higher-order terms (1 — 2y? + 

..). It is clear that no successive approximations can introduce a first- 
order term into 7,,‘. We must conclude, therefore, that the inclusion of the 
heat-capacity effect into our solution for the maximum temperature of the 
micro-meteorite reduces by a first-order term the velocity at which the 
maximum occurs, and correspondingly lowers the atmospheric height and 
increases the density; but it does not change the maximum temperature by 
a first-order term. 

The value of the correction term, y, expressed in terms of the original 
physical quantities involving m/B from equation (18a*), may be written 


y = — (20) 


The quantity y is a maximum at the minimum value of V’.., 11.2 km./sec. 
We may adopt an extremely high value of C, = 0.3 cal./g. for irons or stones 
(C, = 1.3 X 10’ ergs g.), a = 0.9, D = 2.2, T, = 1600°K., and neglect 
To. The resulting value of y appears to be a maximum for any micro- 
meteorites traversing the earth’s atmosphere undamaged. This value is 
y = —0.063. Since the correction term to 7,4 — To‘ is 1 — 2y’,..., the 
maximum correction is — 0.008, less than 1%. 

Of critical importance is the fact that the cos Z term cancels out of y in 
equation (20). This cancellation means that the correction term is inde- 
pendent of the angle at which the micro-meteorite strikes the atmosphere. 
Even though relatively massive bodies may be stopped without melting 
when cos Z is quite small, the effect of heat capacity is no more serious than 
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for a much smaller body siriking vertically. As cos Z decreases, the time 

available for equilibrium increases in such a fashion as to cancel the cor- 

rection term. 

The question of the mean free path of air molecules compared to the 
dimensions of the micro-meteorites must be investigated. Since the 
momentarily trapped air molecules leave the micro-meteorite at the thermal 
velocity, v7, of the surface, their density near the surface tends to increase 
above the ambient air density by a factor of the order of V/vp°. The im- 
pinging air molecules must pass through this haze of escaping gas for a 
distance comparable to the diameter of the body. If the mean free path 
of air molecules in the region surrounding the micro-meteorite near the time 
of maximum surface temperature is greater by a factor of V/vr than the 
diameter of the body, we may be reasonably certain that the present theory 
requires no correction for this effect. There is every reason to believe 
that the mean free paths of the air molecules at the high velocities here en- 
‘countered will be greater than those calculated at normal temperature. 
Hence, if our limit is a safe one for air at normal temperature, it will cer- 
tainly be safe under the conditions of the present theory. 

An inspection of table | in the next section shows that the maximum 
diameter of an undamaged micro-meteorite at vertical incidence velocity 
is of the order of 0.01 em. or smaller. The maximum air density at maxi- 
mum surface temperature is of the order of 6 X 10~* g./em.* At normal 
temperature the mean free paths of air molecules at this density is of the 
order of | cm. Hence the mean free path exceeds the diameter of the 
micro-meteorite by a factor of approximately one hundred times. The 
ratio V/vr is only of the order of ten. Hence the effect of mean free path 
may be neglected unless cos Z < 1/10, in which case 2s may be comparable 
to the standard mean free path. In fact the approximation is rather 
good for radiants to about 1° of the horizon (Z = 89°) or even farther, 
because of the effective increase of mean free path with velocity and be- 
cause the effect of molecular encounters changes the theory slowly as the 
mean free path decreases. 

Since p,, and s vary roughly in the same fashion with |’.,, and since s 
varies as sec Z, the problem of mean free path is a function of Z only. For 
radiants below an altitude of about 1° the calculated maximum dimensions 
are somewhat underestimated in the present formulae because the ac- 
commodation coefficient may be reduced by the shielding effect. 

The time-lag in heat transfer from the front to back surface of the micro- 
meteorite may now be investigated. Suppose that the micro-meteorite 
does not rotate. A slight development of the “cold wave’ approximation 
of heat transfer as used by L. R. Ingersoll and O. J. Zobel‘ can be applied 
to thiscase. A temperature differential over a distance L will be reduced to 
1/e of its value in a time rL* (2//*), where /7* is the heat diffusivity of the 
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material. For ordinary stone //* > 0,005 cm.?/sec. For solid metals its 
value is much larger, With //* = 0.005 and L = 0.01 em. (the largest value 
of 2s for cos Z = 1), the time is 0.03 sec. Unless the meteoric material is a 
much poorer conductor than ordinary stone, the lag in heat transfer will 
not be serious for the largest micro-meteorites at vertical incidence. At 
lower angles of incidence the time of travel increases as sec Z, as does the 
dimension of the body. Hence the time lag increases as sec Z. The lag 
begins to introduce appreciable error for Z > 87°. The effect will reduce 
the radiation rate somewhat below its assumed value in the present theory. 
The calculated critical dimension of a non-rotating micro-meteorite will 
be overestimated at the lowest velocities in the range 87° < Z < 89°, 
If the body is rotating rapidly the error introduced by heat lag will usually 
be small, even to Z = 89°. 

The present theory, however, underestimates the critical dimensions of 
micro-meteorites at extremely large zenith angles because the curvature of 
the Earth is neglected. In no case can the curvature of the path com- 
pletely compensate for this effect in view of the fact that the critical velocity 
is 0.716 V., while the velocity in a circular orbit is 0.707 times the mini- 
muin velocity of fall from infinity. We may roughly define the limits set 
by curvature in the following fashion: while the micro-meteorite is falling 
through a vertical distance 1/, in which the air density increases by a fac- 
tor e, the change in height introduced by curvature, /, must not exceed 6/h, 
where dissmall. The horizontal distance, S, then becomes 


S = (tan Z)/b. (21) 
The corresponding change in height, Ah, arising from curvature over the 
horizontal distance, L, is given to sufficient accuracy by 
Ah = 5*/2R, (22) 
where R is the radius of the Earth. 
If we substitute Ah = 6/b in equation (22) and substitute S from equa- 
tion (21) we find 
tan? Z = 2dRS. (23) 
In an atmosphere of constant temperature gradient, the value of b be- 


comes 


+ ¥) 


7\ p 
Hence the limiting value of zenith distance, Z,,4., from equation (28) 
can be shown to be given by 
tan’ } = . 


7\ p 
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For 6 = 1/100, or 1% accuracy, the values of Zmax lie in the range of 
70° to 77° for Ve having minimum and maximum values, respectively, 
applying to solar-system micro-meteorites. For 6 = 1/10, the corre- 
sponding range in Zax. is 84° to 86° (see next section). 

As a consequence of the above several considerations we may neglect 
effects arising from heat capacity, secondary encounter processes and heat 
transfer in the above explicit theory of undamaged micro-meteorites to 
limits of zenith angle at which curvature of the earth becomes appre- 
ciable. The precision, within the limits of the assumptions, is within 
1° for all micro-meteorites striking at Zenith angle less than roughly 
70% and within about 10% for zenith angles between 70 and 84%. The 
theory probably represents the general order of magnitude to a zenith angle 
of about 87°, when effects of earth curvature become serious for all veloci- 
ties, while the incomplete heat transfer may introduce further errors at 
low velocities. 

Since Earth curvature and curvature of path can be allowed for in the 
general theory where numerical integration is used, this theory can be 
trusted to high precision to a zenith angle of approximately 87° for the 
minimum velocity. At velocities exceeding 20 km./sec., the theory should 
be reliable to a zenith angle of about 89° or perhaps a little greater. 

The fundamental equations to be solved numerically are: 


dV ADpV?, 


2m 
and 
aApVs 
(26 


and 
dt = —dh/(V cos z), (13*) 


in case Earth curvature may be neglected. 

Otherwise / and ¢ must be related by more accurate expressions involving 
earth curvature and the effect of gravity introduced in the motion. The 
atmospheric density, p, must be adopted as a function of height, 4, from 
some standard atmosphere. 

Until a more detailed theory is adopted, such as that by Miss Heineman, 
the drag coefficient, D, may be determined by the expression 

Du 2 8 a (27) 
OV Thy 


‘ In the above equations (11*, 26, 13*, 27) only m, 79, ki, and « may be 
taken as strictly constant in the most general case. ‘The numerical inte- 
gration, however, still remains a relatively simple type of operation. 
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The author has met with little success in calculating the degree of sur- 
face deterioration that may occur to a micro-meteorite at higher velocities 

when atmospheric molecules impinge on the surface with energies of sev- 

eral hundred electron volts. Near the lower limit of velocity, 11.2 km./sec., 
the energies are so small that the deterioration should not be appreciable. 
Near the upper limit of solar-system velocities, 72 kim./sec., the atmos- 
pheric molecules must penetrate deeply below the surface, dissociate and 
damage the lattice structure of the solid. Since the total air mass encoun- 
tered before the velocity becomes greatly reduced is comparable to the 
mass of the micro-meteorite, this damage may be considerable. A crude 
estimate places the critical velocity in the neighborhood of 30 km./sec. 
Above this velocity the damage probably becomes appreciable and near 


TABLE 1 


NUMERICAL VALUES FOR MICRO-METEORITES 


J 
NACA v-2 
mass*/ Mass*/ 
SURFACE, SURPACRE, 
Ve, p (max. T), G./CM.? G./cM.? 


D 104 


11.3 2.09 5.6 10-9 41.4 24.3 
15.0 2.07 24x 19.2 11.8 
20.0 2.05 1-0 >< 107" 8.8 5.6 
23.1 2.04 6.5 x 10°" 6.0 3.9 
25.0 2.04 5.2 « 10°" 4.8 3.2 
30.0 2.08 3.0 x 10°" 2.94 1.99 
40.0 2.02 1.3 x 10°" 1.36 0.96 
50.0 2.02 6.5 107" 0.74 0.54 
60.0 2.02 3.7 X 107"! 0.46 0.34 
70.0 201 24 10° 30 0.23 i 


* The tabulated quantity is radius in microns for a spherical body of density 3(s = 
3m/Bps) at vertical incidence. At other angles of incidence the quantities should be 
multiplied by secant zenith distance of the radiant. See equation (18¢*), et seq., for 
dimensions of non-spherical particles. 


70 km./sec. may become quite destructive. At these velocities, however, 
the maximum dimensions of micro-meteorites are the order of only one 
micron (table 1, above), rather small for easy study at the earth’s surface. 

3. Numerical Results.—The atmospheric density, p», at the maximum 
surface temperature of the micro-meteorite is given by equation (20*). 
Although this equation is derived on the assumption of an atmosphere of 
constant temperature, molecular weight and gravity, nevertheless, it gives 


the correct order of magnitude for p,,. The constants, a, accommodation 

coefficient and, 8, emissivity, we have already adopted as nearly unity and, i 
for lack of better information, equal. If we adopt a maximum temperature, 
Tn.” equal to 1600°K., and 7) = 300°K. then equation (20*) becomes 1 
numerically, 
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Pm = 8.1 X 10°/V.3, (28) 


where c. g. S. units are used. 

Numerical values of p,, are given in the third column of table | for various 
values of V... in the first column. 

The value of the drag coefficient, D, is determined from equation (22*) 
with the substitution of |’,, for V. A mean molecular weight of 26.5 and 
the maximum temperature, 7’,,, are used in calculating, D. The resulting 
numerical equation 1s 


D=2+1.0X 10°/V, (29) 


where is expressed in cm./ sec. 

Numerical values of D are given in the second column of table 1. 

The limiting ratio of mass to total surface area, m/ B, for micro-meteorites 
is calculated only from equation (8), in an atmosphere of constant tem- 
perature gradient. Vertical incidence is assumed so that cos Z = 1. The 
remaining constants depend upon the atmospheric conditions. Two solu- 
tions are given, one based upon the Tentative Standard Atmosphere of the 
National Advisory Committee for Aeronautics’ and the second based upon 
atmospheric results from the V-2 rocket firings. 

The data adopted as a close approximation to the NACA atmosphere 
are: h =OQat 83 km. altitude, = 3.32 10-5 g./em.*, 7) = 240°K., = 
26.5, g = 950¢. g.s.,g = 3.05 K 107-4 °C./em. and vy = 3.65 K 107-8 °C./em. 
The temperature gradient is, therefore, positive with height at a rate of 
3.65 °C./km., beginning at 240°K. at an altitude of 83 km. The resulting 
numerical form of equation (8) is 


m/B = 3.27 X 10% DsecZV,~ *™. (30a) 


Numerical results from equation (30a) for cos Z = | are given in the 
fourth column of table 1. 

Typical results for atmospheric density by measures from V-2 rockets are 
given by E. Durand® of the Naval Research Laboratory. Adopted are: 
h = Oat 100 km. altitude, p, = 1.00 K 10~° g./em.’, 7) = 223°, w = 26.5, 
g = 950 c. g.s.,g = 3.05 K °C./em. and = 5.6 K 1075 °C./em. 
Hence the temperature, beginning at 223° at 100 km., increases with height 
ata rate of 5.6°C./km. The numerical form of equation (30) is 


m/B = 2.55 X 10" D sec Z (30b) 


This equation leads to the results given in the fifth column of table 1. 
A comparison of columns four and five shows that the V-2 atmosphere leads 
to smaller values of m/B than the NACA atmosphere in a ratio of about 
2/3, little dependent upon velocity. This difference arises largely from the 
difference in logarithmic density gradient between the two atmospheres 
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at corresponding values of the density. The critical heights range from 
just below 100 km. to roughly 140 km. 

The calculations of the limits of zenith distance at which the curvature 
of the Earth introduces appreciable error, given by equation (25) and pre- 
sented near the end of the previous section, were derived from the V-2 at- 
mosphere. This atmosphere gives less favorable limits than the NACA 
atmosphere. 

Quantities in table 1 are also presented specifically for V.. = 23.1 km./ 
sec., corresponding to the Giacobinid meteor shower. In the Washington, 
D. C. area, cos Z was about 0.45 for this shower at the time of maximum 
activity. Hence for undamaged iron micro-meteorites of density 7.7, 
the maximum spherical radii should have been 5.2 microns for the NACA 
atmosphere and 3.4 microns for the V-2 atmosphere. For stony-irons of 
density 4.0 these values become 10.0 and 6.5 microns, respectively. For 
long cylinders of circular cross-section, the minor diameter is 4/3 the spheri- 
cal radius, giving a predicted range from 4.2 to 13.3 microns, the range de- 
pending upon uncertainties in the atmosphere and in the density of the 
micro-meteorites. 

Coincident with the shower, allowing for time of fall, H. E. Landsberg’ 
found at Mt. Weather, Va., ‘“‘wedge-shaped opaque” magnetic particles of 
dimensions 40 & 5, 40 * 5 and 100 & 4 microns. These he attributed to 
the shower. The agreement of the minor dimension with the present pre- 
dictions is entirely satisfactory. He also found 3 “round” particles of 
diameter 5 to 10 microns, again in agreement with this theory. On the 
other hand, he found a few “round” particles of diameters from 20 to 40 
microns. Although these are too large for the present theory, there is no 
reason to doubt that larger bodies may have been melted and stopped with 
minor surface vaporization, which would certainly tend to smooth the sur- 
face-—in keeping with Landsberg’s observations. Opik® discusses the 
formation of droplets in meteoroids in considerable detail. 

Landsberg estimates roughly a rate of fall of 3'/2 miles per day for the 
particles of diameter 4 to 10 microns, 13 miles per day for those of diam- 
eter 20-30 microns and 100 miles per day for the one of diameter of 60 
microns. The expected times of fall would be quite variable for the smaller 
particles because of atmospheric turbulence. 

Generally, then, the present theory is in good agreement with Lands- 
berg’s observations. The critical dimensions here calculated are of most 
value in predicting a dimension that divides irregular rough-surfaced par- 


ticles from larger ones that have been melted and partially vaporized in 


passing through the atmosphere. 

It is now clear that the dependence of critical particle size on assumed 
atmospheric conditions is too small to make the study of micro-meteorites 
valuable at this time as a measure of atmospheric density gradient. Never- 
theless, the minor dimensions of the largest rough irregular wedge-shaped 
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pieces observed by Landsberg fall closer to the prediction of the V-2 at- 
mosphere than the NACA atmosphere. It is difficult to know how to 
weight the significance of these results, particularly in view of uncertainties 
in the physical characteristics of micro-meteorites. 

Generally the present theory may be most useful in the following re- 
spects: 

(a) It provides theoretical evidence that micro-meteorites can fall un- 
damaged, hence gives impetus to their collection and study, which should 
lead to measures of space densities for small particles. 

(b) It can aid in distinguishing real micro-meteorites from terrestrial 
objects. 

(c) It can aid in identifying micro-meteorites from two nearly coin- 
cident meteor showers. 

(d) It can assist in determining the original velocity of recognized 
micro-meteorites. 

(e) It may help provide evidence concerning the final dimensions of 
particles producing meteors. 

The collection and laboratory study of micro-meteorites from meteoric 
streams can be invaluable in establishing the nature of interplanetary 
material, specifically cometary material. If the writer’s present working 
hypothesis is correct, viz., that comets were formed at great solar distances 
by the condensation of gas at very low temperatures, then micro-meteorites 
from meteor showers may represent the only laboratory samples of typical 
or atypical interstellar solids. Ordinary meteorites appear to arise from a 
different source.’ 

On a purely speculative level is the possibility of detecting evidence for 
past solar-system catastrophies. If a minor or major planetary disrup- 
tion actually occurred in recent astronomical time, as suggested by C. A. 
Bauer" to account for the meteoritic helium contents, evidence for the 
concomitant micro-meteorites may conceivably be found in the chalk beds 
of the Cretaceous Period or in other geological formations. 

Furthermore, there is the possibility that micro-meteorites in deep oceanic 
sediments (as suggested by the Challenger I:xpedition), in glacial snows 
or even in geological formations (such as chalk beds) may provide a record 
of past astronomical activities. 

The writer earnestly hopes that this paper may encourage the collec- 
tion and study of micro-meteorites. Current work on this problem is being 
conducted by D. K. Norris and F. 5. Hogg,'' while J. D. Buddhue'’ has 
recently summarized his and others’ observation data on meteoritic dust in 
general. 

* This paper was written as a part of the investigation for United States Naval Ord- 
nance Contract-10449-05512. 

1 Whipple, F. L., these PROCEEDINGS 36, 687-695, 1950. 
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MULTIMOLECULAR ADSORPTION OF HORSE ITEART 
METMYOGLOBIN ONTO FILMIS OF BARIUM 


By ALBERT A, Fiskt 


GATBS AND CRELLIN LABORATORIES OF CHEMISTRY, CALIFORNIA INSTITUTE 
OF TECHNOLOGY, PASADENA 4, CALIFORNIAT 


Communicated by Linus Pauling, October 31, 1950 


In a previous paper! it was shown that some protein solutions, when 
adsorbed onto films of barium stearate, followed the classical Langmuir 
adsorption isotherm. The thicknesses of these unimolecular layers were 
calculated by using the Langmuir equation, and it was shown that these 
values for human and horse hemoglobin and for bovine serum albumin 
agreed with values in the literature for the thicknesses of these molecules. 

It was felt advisable to test this technique on another protein of size 
and molecular weight different from the above proteins. Myoglobin was 
chosen because its dimensions have been recently determined by x-ray 
diffraction.’ 

Experimental. The technique was adequately described in the previous 
paper. The thicknesses of the metmyoglobin layers adsorbed onto metallic 
slides covered with an optical gauge of barium stearate were determined 
on the ellipsometer, an optical instrument which measures the ellipticity 
of polarized light reflected from a metallic surface. 

The metinyoglobin solutions were prepared from fresh horse heart by 
repeated precipitations from strong phosphate buffers. The Beckman 
spectrophotometer measurement of a carbonmonoxylated solution did not 
show the presence of any contaminating hemoglobin. The electrophoresis 
of a metmyoglobin solution showed that two components were present, 
the smaller component representing 10°, of the total area. An electro- 
phoresis of a carbonmonoxymyoglobin solution showed that the smaller 
component was present in one-half its concentration in the metmyoglobin 


i 
a 
a} 
1 
: 
| 
* 


Vow. 37, 1951 BIOCHEMISTRY: A. A. FISK 


preparation. A careful separation by electrophoresis was made on a 
metmyoglobin solution, but subsequent electrophoreses of the separated 
main component revealed that the same smaller component was still 
present in 10°; concentration. It was, therefore, assumed that the smaller 
component was probably not an impurity, but was metmyoglobin in the 
dimer or higher form. The adsorption experiments were done on this 
fraction at pH 6.8 in phosphate buffer of ionic strength 0.1. 

Results. The results may be seen in figure 1. The thicknesses in 
figure 1 are not equivalent barium stearate thicknesses, as in the previous 
paper, for refractive index corrections due to the relative humidity of the 
laboratory and the assumed method of packing of the adsorbed myoglobin 
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molecules led to small increases in the thicknesses. It was assumed, on 
the basis of x-ray diffraction data? and the preliminary results, that the 
packing model would be best approximated by assuming the closest pack- 
ing of elliptical cylinders lying on their sides, the voids being filled with air. 
The thickness corrections were calculated by using the derivative of the 
Drude equation,* 
dl 2dn 
n(n? — 1) 


equivalent thickness of barium’ stearate, 
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n = refractive index of barium stearate. 


Interpretation of Results.—The results as shown in figure 1 indicate 
multimolecular adsorption of metmyoglobin onto barium stearate. The 
Brunauer, Emmett and Teller equation‘ for multimolecular adsorption is 


d(Co—C) dnb dnb Cy’ 


concentration of metmyoglobin solution, 
calculated thickness, 
constant, the concentration at infinite thickness, assumed to be 
2.90 g./100 ml., 
= constant related to the heat of adsorption, and 
= thickness of a unimolecular layer of metmyoglobin molecules. 


C) against C/C, using the 
data in figure 1, gives a straight line. From the calculated slope and 
intercept, the thickness of a unimolecular layer of metmyoglobin molecules, 
d», Was determined to be 62 A. 

This value is to be compared with the maximum probable value, 57 A., 
for the diameter of the myoglobin molecule and the value, 9 A., for the 
thickness found by Kendrew’ in his x-ray diffraction study. It would seem, 
therefore, that the metmyoglobin molecule has been adsorbed onto the 
barium stearate surface so that a long dimension, not the thickness, is 
perpendicular to the surface. This case is to be contrasted with the 
previous work! on hemoglobins and bovine serum albumin, where the 
short dimension was perpendicular to the surface. However, since myo- 
globin is known to have but one heme and since this heme is presumed to 
be attached to the side of the molecule,’ it is not surprising to find that 
myoglobin may be adsorbed with the heme in a preferred orientation, 
either against the barium stearate surface or away from it. 

The multimolecular adsorption of myoglobin, in contrast to the uni- 
molecular adsorption of the other protein molecules studied, is not entirely 
unexpected in view of the fact that Taylor® found a considerable amount 
of heme-heme interaction in magnetic susceptibility measurements of 
myoglobin solutions. The ability of myoglobin molecules to associate 
was also apparent in the aforementioned preparative electrophoreses. 

From the above experiment, it is apparent that the determination of 
protein dimensions from measurements of thin films of proteins adsorbed 
onto a barium stearate surface is not limited to molecules in the hemo- 
globin-albumin class. However, it is equally apparent that the interpreta- 
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tion of the results may depend on other information concerning the shape 
and dimensions of tae molecule. 

Summary.—The thickness of a unimolecular layer of horse heart met- 
myoglobin molecules adsorbed onto barium stearate was calculated to be 
62 A. by means of the B.E.T. equation for multimolecular adsorption. 


* This work was supported in part by grants from the American Cancer Society and 
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A COMPARISON OF THE CONTENT OF DESOX YRIBOSEN UCLEIC 
ACID (DNA) IN ISOLATED ANIMAL NUCLEI BY 
CYTOCHEMICAL AND CHEMICAL METHODS 


By CEcILIE LEUCHTENBERGER, ROGER VENDRELY* AND 
COLETTE VENDRELY* 


INSTITUTE OF PATHOLOGY, WESTERN RESERVE UNIVERSITY, CLEVELAND, OHIO 
Communicated by C. W. Metz, November 7, 1950 


In 1948, Boivin, Vendrely and Vendrely' ? presented evidence that the 
amount of desoxyribosenucleic acid (DNA) in nuclei of different tissues 
was constant for the same animal and twice that of the sperm. These 
authors suggested that such a constant relationship applies to all diploid 
animal cells and might be an expression of the genetical equipment of the 
cells. 

While some recent work*~* lends support to this idea, other studies stand 
in contradiction. In 1949, Mirsky and Ris® reported analyses on the con- 
tent of DNA in nuclei of various tissues of mammals which showed a 
quite different relationship to the sperm from that reported by Boivin, 
Vendrely and Vendrely. According to Mirsky and Ris, in calf thymus, 
calf lymph nodes, beef kidney and beef liver, the amount of DNA in 
diploid nuclei in relation to the amount of DNA in the sperm varied from 
2.5 to 3:1, in contrast to the 2:1 ratio found for the same nuclei by Boivin, 
Vendrely and Vendrely. 
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In 1950, Pasteels and Lison,’ working on the rat, reported that the 
DNA in nuclei of diploid liver and pancreas cells showed a relationship to 
the sperm which was only 1.3-1.5:1, which is again a significant departure 
from the ratio found by the French workers. 

Considering that Mirsky and Ris found values of DNA for mammalian 
diploid nuclei which are much too high while Pasteels and Lison found 
values which are much too low in relation to the sperm, such deviations, 
if confirmed, would, of course, be a strong argument against the concept 
of Boivin, Vendrely and Vendrely that the amount of DNA in all diploid 
somatic nuclei of animals is exactly twice that of the sperm. 

A reinvestigation of this particular problem therefore seemed essential 
and it was decided to make a determination of the DNA on the same 
preparations of isolated animal nuclei by two independent, quite different 
methods: the chemical method, in which the amount of DNA per nucleus 
is computed from a large number of nuclei subjected to biochemical 
analysis, and the cytochemical method, which allows a direct estimation of 
DNA in a single nucleus through photometric measurements. In the 
present paper we report the results on the amount of DNA in nuclei of 
these two methods when applied to various tissues of beef and rat. 

Material and Methods.—Liver, kidney, spleen and sperms were obtained 
from freshly killed beef and rats, the nuclei isolated and the amount of 
DNA determined chemically as previously described by Vendrely and 
Vendrely.? For the cytochemical analysis, the same suspension of isolated 
nuclei was used for smears, which were fixed in Carnoy and 50% Formalin 
and on which the Feulgen reaction was carried out. The DNA was esti- 
mated by photometric measurements with an apparatus similar to that 
of Pollister and Moses,* and by the method as described by Schrader and 
Leuchtenberger.4. The cytochemical analysis was also done on sections of 
the same organs fixed in Carnoy and 50°) Formalin. 

Results.-In the following table the results of the cytochemical and 
chemical analysis of DNA in isolated nuclei of different tissues of beef and 
rat are presented. It should be pointed out that the amount of DNA per 
nucleus found-by the cytochemical analysis is a mean value expressed 
in arbitrary units (as previously deseribed),‘ while the amount of DNA 
per nucleus by the chemical analysis is given in absolute amounts in 10~® y, 

If we consider first the cytochemical analysis of the DNA of different 
tissues of the beef, it is evident that the liver, spleen and kidney contain 
approximately the same amount of DNA per nucleus and twice that of the 
sperm. The proportion is also maintained if we look at the data obtained 
by the chemical analysis on the same suspension of nuclei. If we now 
examine the cytochemical analysis of the DNA in the tissues of the rat, it is 
evident that in the liver the nuclei fall into three classes with values of 
DNA per nucleus of 3.5, 6.6 and 12.1 (arbitrary units), respectively. 
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Furthermore, these values represent nearly exact multiples of each other. 
The picture is quite different if we consider the amount of DNA of the 
nuclei of the kidney of the rat. Here only a mean value of 3.3 per nucleus 
was observed, which corresponds to the lowest class of DNA found in liver 
nuclei. 

In the chemical analysis, all the size classes of nuclei in the liver are, 
of course, thrown together and the results of such lumping is expressed by 
the fact that the amount of DNA per liver nucleus is considerably higher 
than that of the kidney. 

Discussion of the Results.—In evaluating the data it should be pointed 
out that the cytochemical values given in the table for the amount of 
DNA per nucleus of a tissue actually represent a mean of a number of 
individual DNA values. While this mean amount is constant and the 


TABLE 1 


COMPARISON CF THE AMOUNT OF DESOXYRIBOSENUCLEIC AcID (IDNA) IN ISOLATED 
NuCLEI OF DIFFERENT TiSSUES OF BEEF AND Rat By CHEMICAL AND CYTOCHEMICAL 


ANALYSIS 
~—-CYTOCHEMICAL ANALYSIS 

MEAN AMOUNT OF CHEMICAL ANALYSIS, 
TYPHK OF TYPE OF NO. OF NUCLEI DNA PER NUCLEUS AMOUNT OF DNA PER 
ANIMAL TISSUE MEASURED (ARBITRARY UNITS) nucLEus (10 
Beef Liver 30 3.3 = 0.11 6.8 
Beef Spleen 23 3.0 * 0.08 7.0 
Beet Kidney 26 3.3 = 0.10 6.4 
Beef Sperm 25 1.6 + 0.03 3. 3° 
Rat Liver 15 3.3 = 0.08 
Rat Liver 21 6.6 + 0.24 8.2 
Rat Liver 7 12.1 = 0.50 
Rat Kidney 38 3.3 = 0.06 5.9 


® This value is taken from previous analyses by Vendrely and Vendrely.* 
Nore: The amount of DNA determined by cytochemical analysis is given in arbitrary 
units. 


same for diploid nuclei of different tissues for the same aninial there is a 
variation from nucleus to nucleus within each tissue, sometimes as high 
as 50%. For example, the beef spleen which gave an average amount of 
DNA of 3 per nucleus, similar to the average amount for beef kidney and 
beef liver, showed individual nuclei with values of DNA ranging from 2.7 
to 4.0. Whether such differences are due to errors inherent in the cyto- 
chemical photometric method as previously discussed,® or are due to 
biological variations in the amount of DNA among individual nuclei of a 
tissue cannot be decided at the present time. 

If we compare now the results obtained by the cytochemical and the 
chemical analyses of DNA in the nuclei of different tissues of the beef, 
it is evident that by both analyses the diploid nuclet of the somatic tissues 
show twice the amount of DNA of the sperm. Therefore, the unavoidable 
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conclusion follows that constancy of the DNA in nuclei of the beef, as 
first reported by Boivin, Vendrely and Vendrely in 1948, holds true, even 
when individual nuclei are analyzed for their content of DNA (as done by 
the cytochemical studies), instead of analyzing a large mass of nuclei and 
computing the value for one nucleus. These results, therefore, do not 
confirm the measurements of DNA of mammalian nuclei as reported by 
Mirsky and Ris. “According to the last-named authors, the higher values 
of DNA in the beef liver (which represent three times the value of the 
sperm—-namely, 8.4 to 2.82) cannot be explained by the presence of poly- 
ploid nuclei, for they state: ‘‘Polyploidy is not the explanation of the high 
values found for beef liver, for in counting liver nuclei the double sized 
nuclei were counted separately and the necessary correction was made.” 

In view of the fact that the results which we obtained in the present 
study by both eytochemical and chemical analyses are in complete agree- 
ment with the findings previously obtained on beef by Vendrely and 
Vendrely, the discordant findings of Mirsky and Ris remain unexplained. 
_ It might be of interest to state that in the cytochemical analysis of DNA 
of the nuclei of the beef liver only two nuclei among thirty-two examined 
showed amounts of DNA which were twice that of the values reported. 
A similar relationship holds true for spleen and kidney. It should be pointed 
out, however, that occasionally nuclei were found which had a considerably 
lower value than 3.3, especially if nuclei of very small size with a pycnotic 
appearance were measured. Whether these are degenerating, dying nuclei 
so familiar to every pathologist, or whether the lower value is due to errors 
in the photometric measurements cannot be decided at this moment with 
certainty. We are inclined to believe that even in the normal organs there 
are always a small number of degenerating nuclei present which gradually 
lose their DNA as demonstrated for pycnotic nuclei by Leuchtenberger.® 

Turning now to data on rat nuclei, the finding of three distinct classes 
with multiple values of DNA, 3.3, 6.6 and 12.1 respectively in the analysis 
of isolated liver nuclei is in good agreement with the three classes of DNA 
values found by Ris and Mirsky” in sectioned material of rat liver. A 
similar agreement between amounts of DNA found by the cytochemical 
analysis in isolated nuclei and tissue sections has been found also in our own 
studies. 

Pasteels and Lison’ also report three classes of DNA values in the rat 
liver nuclei in sections but each of these classes has a value which is about 
30% too low, if compared with the results in other organs. While kidney, © 
for instance, has one class of DNA with a value of 227, the liver has the 
three size classes—-156, 347, 772. Hence the measurements in the kidney 
find no equivalent in the liver nuclei where the minimum value is signifi- 
cantly lower, while the two higher values stand in no multiple relation 
to the single class represented by the kidney nuclei. 
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We cannot confirm these results because, although we found also only 
one class of DNA in kidney in contrast to the three classes in the liver as 
did Pasteels and Lison, the amount of DNA in the nuclei of the kidney of 
the rat is approximately the same as the lowest amount of DNA in the 
liver nuclei. 

If we now compare the results obtained by our cytochemical analysis 
of DNA in nuclei of rat liver and rat kidney with the chemical analysis on 
the same material, it is evident that the higher value of DNA in the liver 
nuclei—8.2 X 10~*y as compared with the lower value of 5.5 & 10~*y for 
the nuclei of the kidney—is an expression of the polyploidy in the liver. 
We feel that these studies demonstrate rather well the importance of the 
simultaneous application of the cytochemical and chemical method for the 
analysis of DNA on the same material. Only by the cytochemical method 
was it possible to recognize in the rat liver three different types of nuclei 
which were characterized by three different amounts of DNA (multiples 
of each other). In contrast the result of the chemical method represents 
an average value of three types of nuclei with different amounts of DNA 
and is, therefore, not the true quantity of DNA of the individual nuclei in 
the liver. This does not by any means lessen the importance of the 
chemical analysis of DNA in diploid nuclei, particularly if we keep in mind 
that if the chemist presents a value of DNA for a nucleus it is based on 
much larger samples of nuclei (billions), while the cytochemist—even if 
he were able to study 1000 nuclei of one tissue—always analyzes only at 
comparatively small population of nuclei in respect to the whole tissue. 
There seems to be no doubt as to the desirability of employing both methods 
simultaneously on the same material, especially in cases involving poly- 
ploidy or in tissues which show pathological changes. 

Summary and Conclusions.—DNA was determined in the same prepara- 
tion of isolated animal nuclei by a cytochemical and a chemical method. 
Contrary to the findings of Mirsky and Ris, we found by both methods that 
diploid nuclei of various tissues of beef contain the same amount of DNA 
and twice that of the sperm. The constancy in the amount of DNA was 
also observed in the diploid class of liver and kidney nuclei of the rat in 
contrast to the findings of Pasteels and Lison. 

This analysis, based on two different and independent methods, thus 
gives firm support to the concept of Boivin, Vendrely and Vendrely, that 
mammalian diploid nuclei of various tissues contain a constant and charac- 
teristic amount of DNA and that this amount is double that which is 
carried by the sperm."! 


* Faculté de Medecine, Institut de Bacteriologie, Strasbourg, France. 
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STRUCTURES IN Ek. COLI RESEMBLING CIITIROMONEMATA 
By ALFRED MARSHAK 


New York UNtiversity, COLLEGE OF MEDICINE AND THE DIVISION OF TUBERCULOSIS, 
U. S. HEALTH SERVICE 


Communicated by Barbara McClintock, November 11, 1950 


Although bodies having the characteristics of nuclei have been described 
in studies of E. coli with the light microscope and with the electron micro- 
scope, the internal structure of the nuclei or of the chromosomes has not 


been observed hitherto.'~‘ The observations reported here show the 
presence of structures resembling chromonemata, and chemical analyses of 
the particulate fraction in which the apparent nuclear bodies are concen- 
trated fulfill expectations from previous studies.* 4 

Bacteria in the early log phase of growth were suspended in 5% citric 
acid and immediately subjected to sonic vibration (9 ke.)® for 45 minutes 
and then subjected to differential centrifugation. A heavy black pig- 
mented fraction was obtained at 900 g. and a light unpigmented one at 
3600 g. Very little additional solid material was collected at 55,000 g.® 
The average size of the heavy particles was 0.96 ¢ 0.002 u X 0.63 10.002 u, 
which was approximately the size of the Feulgen-positive bodies of the 
unbroken cells. 

The heavy fraction was suspended in distilled water and a drop of the 
suspension on a collodion film supported by wire mesh was frozen on a block 
of solid CO,. The frozen material was then evacuated and shadowed with 
chromium. When examined in the electron microscope many of the 
particles showed rather regular peripheral lobulations, suggesting the 
presence of a banded or helical structure. In many instances, helical 
structures whose coils corresponded with the peripheral lobulations could 
be observed (figure 1). Most of the particles observed gave some indica- 
tion of the presence of helices although these structures appeared much 
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more clearly in some than in others. In some particles, a system of at least 
two helices could be seen and the constituent strands of each helix were 
observed to be themselves helically coiled (figure 2). Occasional particles 
were found whose widths were approximately one-half those of the larger 


be 


FIGURE 
Particles from the heavy fraction showing peripheral lobulation and coils 
of helices. Chromium 17,500 . 
FIGURE 2 
Helical coiling of individual strands composing the larger helices. Chro- 
mium 52,500 
FIGURE 3 
Particle of one-half the usual width showing the helix composed of 2 strands 
each of which is itself helically coiled. Chromium 52,000 *. 


particles (figure 3). They could be seen to contain a rather regular helix, 
the gyres of which were made up of two strands. Where the individual 
strands could be cleariy distinguished, they could also be seen to be helically 
coiled, with a pitch approximately |; that of the larger helix. 
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After extraetion with cold trichloracetic acid and with fat solvents the in- 

tact bacteria and the particulate fractions were digested with perchloric 

acid and the purine and pyrimidine content determined by chromatography 
and spectrophotometry.” > The results obtained are given below: 


Nucrerc Acip Bases (uG./McG. Dry MATTER) 


| ADENINE GUANINE CYTOSINE URACIL THYMINE 


Intact bacteria 18.6 24.3 13.9 7.39 3.32 
Heavy particles 10.6 15. 12.1 4.46 2.34 
Light particles 19.9 25 14.8 8.76 .85 


Each value in the table represents a single determination and is subject 
to an error of approximately 10°). Absorption spectra were taken of each 
recovered base, and comparison with the curves for the pure bases showed 
good agreement except in the case of thymine. The curve for thymine from 
the intact bacteria gave a fair approximation to the curve for pure thymine, 
but the curves for material recovered from the particulate fractions devi- 
ated markedly. 

The contamination has the effeet of increasing the apparent concentra- 
tion of thymine but the extent of contamination cannot be determined with 
accuracy. The thymine values must therefore be considered to be unreli- 
able. The reasons for the difference in chromatographic behavior of the 
digests from the whole bacteria and the particles are unknown, 

The concentration of each base other than thymine is lower in the heavy 
particle fraction than in either the light particles or the intact bacteria. 
The sum of the bases is lower in the heavy particle fraction than in either 
the intact bacteria or the light particle fraction by 34-350). 

Discussion... Previous investigations have shown that the cytoplasm 
of E. coli is more dense than the nucleus,* and that the material responsible 
for the greater density can be removed by procedures for extracting ribo- 
nuclere acid.! If it is postulated that the heavy particle fraction represents 
a concentration of nuclear material as suggested by the dimensions of the 
particles and their internal morphology, a lower concentration of nitroge- 
nous base would be expected. This expectation was fulfilled. A greater 
proportion of thymine to the total base content would also be expected. 
However as indicated previously the thymine data are unreliable. As far 
as they go the chemical data support the hypothesis that there is a concen- 
tration of nuclear material m the heavy particle fraction. 

The chromatinic bodies described by Hillier, Mudd and Smith® are fre- 
quently bipartite or bifureate. Rough measurements from their illustra- 
tions give an approximate value of 0.3 « for the width of a single limb of these 
bodies. This may be compared with the width 0.29 w obtained from 
measurement of particles containing a single large helix (figure 3) and the 
width 0.65 yu for particles containing two such helices (figures 1 and 2). 
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If, as has been previously suggested, each lirab of the bipartite body be 
assumed to be a chromosome, then the above calculations suggest that the 
helical structures may be chromonemata. The bipartite nature of the 
helix and the presence of ‘‘secondary coiling’’ of the constituent helices are 
exactly analogous to conditions found in chromosomes of higher plants.® '° 
On the basis of both the morphological and the chemical evidence therefore 
it is concluded that the coiled structures here described represent chromon- 
emata of the chromosomes of E. coli. 

Coiled structures have been found in biological materials which are not 
related to chromosomes. ‘The double helix of elastic tissue differs morpho- 
logically from the E. coli structures in that the pair of filaments is not itself 
helically coiled, and also chemically since the coiled fibers of elastic tissue 
disintegrate at a pH lower than 5.7.'' The helices of the sperm neck-piece 
and middle-piece show no “secondary coiling,” the latter occurring in 
groups of 3. Both must be impregnated with chromic or osmic acids to be 
detected, suggesting a high lipid content." 

Summary.—E. coli was fragmented by sonic vibration in the presence of 
citric acid and the fragments separated by centrifugation into a light and a 
heavy fraction. Analysis of purine and pyrimidine content showed the 
latter to represent a concentration of nuclear material. The dimensions of 
the fragments were comparable to those of the chromatinic bodies of the 
same cells. Particles in this fraction showed dual helically coiled structures 
resembling chromonemata of higher plants and differing in morphology 
from coiled structures of non-chromosomal materials previously described. 
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PARTIAL DIFFERENTIAL EQUATIONS AND GENERALIZED 
ANALYTIC FUNCTIONS. SECOND NOTE 


By LipMAN BERS 
SyRACUSE UNIVERSITY AND THE INSTITUTE FOR ADVANCED STUDY 


Communicated by Marston Morse, November 4, 1950 


1. Introduction.—This note (an abstract of results to be published in 
full elsewhere) contains a reformulation and extension of the theory 
outlined in the preceding note under the same title! (henceforth referred 
toas 1). The new theory is consequently based on a complex differentia- 
tion process (and on the concomitant integration) and is therefore a true 
generalization of classical function theory. It applies to all elliptic linear 
homogeneous partial differential equations of second order (and to systems 
of first order equations) for unknown functions of two variables. The 
new theory contains as special cases and unifies the results of I on self- 
adjoint equations, the earlier more special theory of sigma-monogenic 
functions by A. Gelbart and the author,’ and the very ingenious extension 
of this theory proposed by Markushevitch.* 

Notations: Eo denotes the finite z-plane: |z| = |x + iy| < +, E 
the whole function theoretical plane (Riemann sphere): z| s+. 
D the closure of a domain D. Functions of the real variables x and y are 
written as functions of z = x + zy, without implying analyticity. If fisa 
complex number, / denotes its complex conjugate. Partial derivatives are 
denoted by subscripts; in particular, 2w, = w, — tw,, 2w; = w, + 1wy,. 
A (complex-valued) function w(z) belongs to class @ (class 3C) if it is con- 
tinuous (satisfies a Hélder condition). If w has partial derivatives with 
respect to x and y of order n, which are of class @ (of class 3C), then w is of 
class @” (of class ”). w(z) is said to have property @ at z = © if w(1/z) 
has property P at z = 0. 

2. Differentiation and Integration. Pseudo-analytic Functions.--Two 
functions, F(z) and G(z), of class %' defined in a domain Dy ¢ Ep are said 
to form a generating pair in Dy if Im (FG) > 0. Our generalization of 
function theory is based on assigning to a generating pair (F, G) the part 
played in the classical theory by (1, 2). 

A function w(z) defined in a domain D, D ¢ Do, admits the unique 
representation 


w(z) = F(z) + ¥(z)G(z), (1) 


where ¢g and y are real-valued. If for some 2 € D the limit 


w(2o) = lim w(z) — F(z) — p(20)G(2) 
Zo 
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exists and is finite, w(z) is said to possess at 2 the (F, G)-derivative w = 
dip.gqw/dz. It turns out the existence of w in a domain has as far-reaching 
consequences as that of the classical derivative. 

A function possessing an (F, G)-derivative at all points of D is called 
regular (F, G)-pseudo-analytic of the first kind (or, more briefly, pseudo- 
analytic) in D. If (1) is such a function, we call w(z) = g(s) + i(z) 
pseudo-analytic of the second kind. The correspondence between w and 
w is one-to-one, and the two functions should be thought of as two repre- 
sentations of the same mathematical entity. The class of pseudo-analytic 
functions is closed under addition, multiplication by real constants, and 
uniform convergence. 

TuHeEorEM. [f w(z) 1s pseudo-analytic of the first kind, w(z) and w(2) are 
of class H'. If w(z) is pseudo-analytic of the second kind, w(z) 1s of class 
K?, and the mapping w = w(z) is an interior, sense-preserving, quasi-con- 
formal transformation. 

An isolated singularity 2 of a single-valued pseudo-analytic function 
w(z) is called a pole if w(z) ~ a(z — 2)", 2 — 2, wherea # Oand n isa 
positive integer, an essential singularity if w(z) comes arbitrarily close to 
every complex number as 2 —> 2p. 

THEOREM. An isolated singularity of a single-valued pseudo-analytic 
function is either essential, or removable, or a pole. 

Let W(z) be a continuous function, and I a rectifiable curve leading from 
2 to 2. The (F, G)-integral of W over I is defined by the equation: 


GW dz : FW dz 
= F(z)Re Gere f == (8) 


W is called (F, G)-integrable in a domain D if its (F, G)-integral over a closed 
curve I vanishes whenever I’ is homologous to zero in D. 

THEOREM. A continuous function is (F, G)-integrable if and only if tt 1s 
an (F, G)-derivative; 


= — F(21) — p(20)G(2). (4) 


A generating pair (F\, G;) is called a successor of (F, G) if (F, G)-deriva- 
tives are (Fi, G,)-pseudo-analytic, and (F\, G,)-pseudo-analytic functions 
are (F, G)-integrable. (F, G) is then called a predecessor of (Fi, Gi). 

Tueorem. Let (F, G) be a generating pair in Do. In every bounded 
domain D,, D; € Do, (F, G) possesses a successor and a predecessor. 

A sequence of generating pairs {(F,, G,)}, v = 0, 1, #2, ..., is called 
a generating sequence if (F, , ;,G, + ,) isa successor of (F,,G,). A generating 
pair (F, G) in Do can be embedded (in every bounded domain D,, D,; ¢ Dy) 
in a generating sequence |(F,, G,)} such that (/, Go) = (F,G). With 
respect to such a sequence an (F, G)-pseudo-analytic function w(z) has 


= 
/ 
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derivatives of all orders: 
=w, w"t” = n= 1,2,... (5) 
The sequence {w"(z)} determines the function w(z) uniquely. 

Two generating pairs, (F, G) and (F, G) are called equivalent if F =. 
anuF + ayG,G = anF + anG, the aj being real constants. A generating 
sequence {(F,, G,)} is called periodic if there exists an integer » > O such 
that (F, 4,, G,+,) is equivalent to (F,, G,). In a bounded domain every 
generating pair can be embedded in a non-periodic generating sequence. 
I conjecture, but have not proved, that generating pairs cannot be em- 
bedded in periodic sequences, except in special cases. 

Remark: Since (F,G)-differentiability is a conformally invariant 
property, pseudo-analyticity can be defined on a Riemann surface. 

3. Formal powers.—A generating pair (F, G) in Ep is called a generating 
pair in E if F(z) and G(z) are of class 3C' at z = ©, and Im{F(@)G()} >0. 

Turorem. A generating pair (F, G) in E can be embedded in a generating 
sequence {(F,, G,)} in Ey in such a way that the functions F,(z), G,(z) are 
continuous atz = ©. The generating pairs (F,, G,) are determined uniquely, 
up to equivalences. 

In this and the next section we consider a fixed generating pair (F, G) 
in E embedded in a sequence described in the theorem. This involves no 
serious loss of generality, since given a generating pair in Dy ¢ Ep, we can 
find a generating pair in E which coincides with the given one on a pre- 
assigned closed subdomain D, ¢ Do. 

Let a and 29 be complex numbers, a # 0, and ra real rational number. 
If r # 0, set |r| = p/q, p and q being relatively prime positive integers; 
ifr = 0, set p = 0,q = 1. An (F,, G,)-pseudo-analytic function w(z), 
q-valued and regular for 0 < |z — 20| < + and such that w(z) ~ a(z — 
zo)’, 2 > &, and w(z) = O({2|"), 2 > ©, is called a (global) formal power 
and is denoted by 7") (a, zo; z). The corresponding function of the second 
kind is denoted by +7” (a, 20; 2). Fora = 0, we set Z" (0, 20; 2) = 0. 

Turorem. The formal powers Z"? (a, 20; 2) exist and are unique (for 
allv,r,a,%). Ifa # 0, then the mapping w(z) = *Z'? (a, 2, 2) is a homeo- 
morphism of the q-times covered 2-plane (with branch-points at 29 and @) onto 
the p-times covered w-plane (with branch-points at 0 and ~), Z'? (a, 20, 2) 
is a continuous function of 2. The formal powers satisfy the relations; 


Z (da + ub, 20; 2) = AZ! (a, 20; 2) + (a, 20; 2) 
where \ and yw are real constants, and 
(a, %; 2) = a, Z0(a, %; 2)/dz = 1} (a, 2). (7) 


It follows from (7) that once the sequence {(F,,G,)| is known, the formal 
powers Z™,n = 1, 2, ..., can be obtained by successive integrations. 


| 
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The same procedure can be carried out starting with an arbitrary generat- 
ing sequence defined in a bounded domain Do. In this way we obtain 
local formal powers, which have some properties in common with the global 
powers. (Local formal powers with negative and/or fractional exponents 
can also be defined, but not uniquely.) 

.4. Applications of Formal Powers.—-With the aid of the global formal 
powers we can reestablish all results from I §§5-15 in the present more 
general circumstances. Thus we obtain Taylor and Laurent expansions 
for single-valued and multiple-valued functions, interpolation theorems 
for rational (regular but for a finite number of poles) pseudo-analytic 
functions, the analogs of the fundamental theorem of algebra, of Cauchy's 
integral formula,‘ of Weierstrass’ convergence theorem, of Runge’s 
approximation theorem, etc. We also state a result which describes the 
totality of single-valued pseudo-analytic functions. 

Let w(z) be single-valued and (F, G)-pseudo-analytic in a domain D, 
and let f(z) be single-valued and analytic in this domain. We say that 
w(z) and f(z) are similar if there exists a continuous function x(z) in D, 
such that w(z) = e* f(z). 

THEOREM. Every single-valued analytic function f(z), z¢€D, has a similar 
(F, G)-pseudo-analytic function w(z), and vice versa. If one of the functions 
is given, the other may be chosen so that {[w(z)/f(z)] —> 1 as z —> Zo, 2 being a 
prescribed point of D. 

This theorem implies, among other things, the existence of pseudo- 
analytic functions with infinitely many prescribed isolated zeros and poles. 

5. Pseudo-analytic Functions and Partial Differential Equations.—In 
order to express pseudo-analyticity by differential equations we associate 
with every generating pair (/’, G) the functions 


a = —(FG, — F,G)/(FG — FG), b = (FG; — F.G)/(FG — FG), (8) 
A = —(FG,~ FG)/(FG — FG), B= (FG, — FG)/(FG — FG), (9) 
a — i8 = 4(F/F)(FG, — F,G)/(FG — FG), (10) 


a and being real-valued. 

Lemma. Let Dy and D, be domains such that D, ¢ Dy ¢ Ep. (a) Given 
two functions of class #, a(z) and b(z), defined in Do, there exists a generating 
pair (F, G) in D, such that (8) holds. (b) Given two real-valued functions of 
class #, a(z) and B(z), defined in Do, there exists a generating pair (F, G) 
in D, such that (10) holds. 

The statements (a) and (b) remain valid for D = Dy = Ep if a, b, a, 8 = 
0 

TuHroreM. A function w(z) = F(z) + ¥(2)G(z) of class is (F, G)- 
pseudo-analytic if and only if 
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w; = aw + bw (11) 
1.e., of and only if 
Fe, + Gy = 0. (12) 
If this condition is satisfied, then 
w= w,— Aw — Bw= Fo, + Gy, (13) 
and 
(w), = aw - Bw. (14) 


(Note that (14) and part (a) of the lemma imply the existence of 
successors. ) 
Eliminating from (12) the function g we obtain the equation 


Wrz + Vy + Op, + By, = 0 (15) 


where a and @ are given by (10). 

Example: Let o(z) be a positive function of class K', and set (F, G) = 
io"), Then the (F, G)-pseudo-analytic functions of the first 
(second) kind are the reduced pseudo-analytic (pseudo-analytic) functions 
of I. If o = o;(x)o2(y), then the functions of the second kind are sigma- 
monogenic.’ In this case a generating sequence can be obtained by setting 
(F2, Ga) = (F,G), (Fa +1 Ga 41) = io"), where = o2(y)/o1(x). 


Consider now the general linear homogeneous partial differential equa- 
tion (with real coefficients) : 


Aer + ry + Andy + Br + By, + TY = 0. (16) 


We assume that the coefficients A,, are of class 1', the coefficients B,, of 
class i, and that the ellipticity condition (A;,;A% — A ip > 0) is satisfied in 
the domain considered.’ By introducing new independent variables (16) 
can be reduced to the form 


Vir + Vy + abe + Bh, + W = 0. (17) 


If y = 0, this equation is at the form (15) and it follows from statement 
(b) of the lemma that, for an appropriate choice of (7, G), the solutions 
of our equation will be the imaginary parts of (F, G)-pseudo-analytic 
functions of the second kind. If y # 0, we consider equation (17) in a 
domain where it possesses a positive solution Yo, and introducing the new 
unknown function y = ~/po reduce the equation to the form (15). 

Consider next the general system of linear homogeneous partial differen- 
tial equations (with real coefficients). 


+ + Bats + Bab, + Tne + = 0,7 = 1, 2. (18) 
This system is elliptic if the roots of the equation det (A, — ABy) = 0 
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are complex. We assume that this condition is satisfied in the domain 
considered, and that the coefficients A,, By are of class i', the coefficients 
l', of class H. By introducing new independent variables (18) can be 
reduced to the form 


= He + ofy + + Gy = + Hy + + vow (19) 


where o has constant sign. Without loss of generality we assume that 
o>0. Set y? = yi + yatountye-. Iy=0, equations (19) are 
equivalent to equation (12) for F = 1,G = — ++ 10. Solutions (¢, p) are 
then real and imaginary parts of (F, G)-pseudo-analytic functions of the 
second kind. If y # 0, we reduce (19) to a system of the same form with 
y = 0, by introducing new unknown functions ¢, w such that g=¢eot 
ow, = + poy, where (¢, (¢2, are fixed solutions of (19) with 
¢giv2 — gyi > O. The existence of such solutions follows easily from 
statement (a) of the lemma. 

On the other hand, we can introduce new unknown functions u = g — ty, 
v = oy. Then w = u + wW satisfies an equation of the form (11). It 
follows from (a) that, for an appropriate choice of (F, G), solutions (x, v) 
are real and imaginary parts of (F, G)-pseudo-analytic functions of the 
first kind. 

Thus the theory of pseudo-analytic functions bears the same relationship 
to the general theory of elliptic equations as the theory of analytic functions 


does to that of the Laplace equation. 


1 Bers, L., Proc. Natl. Acad. Sci., 36, 130-136 (1950). 

2 Bers, L., and Gelbart, A., Quart. Appl. Math., 1, 168-188 (1943); Trans. Am. Math. 
Soc., 56, 67-93 (1944); Ann. Math., 48, 342-357 (1947); Bers, L., Am. J. Math., 
72, 705-712 (1950). Fora generalization sce Lukomskaya, M. A., Doklady Akad. Nauk, 
73, 895-888 (1950). 

3 Reported by I. G. Petrovskii, Uspekhi Mat. Nauk, 1, 44-70 (1946). 

‘ For other generalizations of Cauchy's formula see Bergman, S., Trans. Am. Math. 
Soc., 62, 452-497 (1947); Polodii, G. N., Mat. Shornik, 24 (66), 375-384 (1949); Sabat, 
B. V., Doklady Akad. Nauk, 69, 305-308 (1949). 

5 We do not consider here the interesting case of equations possessing parabolic lines; 
cf. reference 2, and Weinstein, A., Trans. Am. Math. Soc., 63, 342-354 (1948); Ouart. 
Appl. Math., 5, 429-444 (1948). 
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THE HARMONIC OPERATOR FOR EXTERIOR DIFFERENTIAL 
FORMS* 


By MatTruew P. GAFFNEY 
NORTHWESTERN UNIVERSITY 
Communicated by M. H. Stone, October 30, 1950 


Harmonic integrals were introduced by Hodge.’ Kodaira® * and 
Bidal and de Rham* ° independently introduced the harmonic operator 
A(= dé + 4d); one of their principal uses for A was to study harmonic 
integrals. A is defined for exterior differential forms on a Riemannian 
manifold of class C’; for functions it reduces to the Beltrami operator. If 
we properly select its domain, A can be regarded as a linear transformation 
on a Hilbert space——the Hilbert space is obtained by suitably completing 
the class of continuous square integrable forms—and the first object of 
our work is to develop the properties of A from this point of view (thus 
avoiding the use of integral equation theory). Our second—and _ prin- 
cipal——object is to extend our knowledge of A in the non-compact case 
(in which the integral equation theory does not work well). Here the 
most important question is whether or not the closure, A, is self-adjoint; 
for only when it is can the powerful spectral theorem be applied.® 

The domain of d (respectively 6) we take to be the set of all @ in C! 
such that both a and da (respectively da) are square integrable. Similarly, 
the domain of A is the set of those a such that a, da, and 6a are all in the 
domains of d and 6. With this definition A will not, in general, be sym- 
metric; we are therefore led to restricting the class of manifolds considered. 
We demand that d and 6 (with the domains provided) be adjoint [(da, 8) = 
(a, 68)]. When this condition is satisfied we say that the manifold has 
“negligible boundary.” A is then easily seen to be symmetric; our main 
result is that for manifolds with negligible boundary A is self-adjoint. 

If a manifold is compact it has negligible boundary. Another sufficient 
condition is that there exist a function s with the following properties 
(which are suggested by a consideration of the function equal to the 
distance from a variable point to the origin in euclidean space): 

I. Outside some fixed compact set K the function s is of class C'. 

II. There exists a number & such that for all t 2 & the set S, = [P: 
s(P) S t] is a compact set containing K. 

III. Outside K the gradient of s is never zero and the magnitude of the 
gradient is bounded. 

These properties can be weakened. We can also give a different suffi- 
cient condition which applies, for example, to the euclidean plane with a 
point removed. 

The first step in proving self-adjointness is to prove that d = 5*: for 


— 
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every a in the domain of 6* we must find a sequence <a, > of forms in the 
domain of d such that in the Ly norm [norm @ = (a, a)} a, approaches 
a and da, approaches 6*a. We first localize the problem to a single 
coordinate neighborhood by the use of a partition of unity, | = 2 @,;. 
To use the partition we must prove that if a is in the domain of 6* so also 
is da. (While this follows readily, the corresponding fact for A* offers 
considerable difliculty.) To produce the approximating sequence <a», > 
we use the Friedrichs mollifier;”: * this is an integral operator on E” whose 
kernel is an approach to the Dirac delta function. On the real line (for 
example) the kernel is defined by j,(x — y) = (== 2), where f(t) is a 
C” non-negative function which is zero outside the interval [—1, 1] and 
whose integral over this interval is 1. For a form a@ (which by localization 
vanishes outside a coordinate neighborhood) we first fix the coordinate 
system and then apply the mollifier to the individual coefficients. Denot- 
ing this operator by J/,, we prove that if ais an L, form then J, approaches 
a. Furthermore, if @ is also in the domain of 6*, then d/,a = J,6*a and 
consequently approaches 6*a; this proves that d = 6*, and it follows that 
6 = d*. By a standard theorem dé and 6d are therefore self-adjoint, and 
with the aid of the relations dd = 66 = 0 it can be shown that dé + 6d 
is self-adjoint. 

We must next show that A = dé + 6d. First we prove that if @ is in 
the domain of d* then 6/,a approaches d*a; since 6 and J, do not commute 
this requires a special argument, but the methods of Friedrichs provide 
the necessary tools. Then if @ is in the domain of 6d, we have 6d/,a@ = 


- 6J.da approaches 6da—note the important role played by the commuta- 


tivity of d and J,. Next it is shown that if 8 is in the domain of d6 + bd 
so also is $,8; this proof uses the facts concerning 6/, and also the in- 
equality we give below. Finally, techniques are developed to prove that 
AJ,(¢:8) approaches (dé + 

The inequality mentioned is of independent interest. For a particular 
coordinate neighborhood denote by D(a) the integral of the sum of the 
squares of all possible first partial derivatives of the coefficients of a. 
Then for every point P there is a coordinate neighborhood A such that 
for all C? forms vanishing outside A 


D(a) = C{(da, da) + (6a, ba) + (a, a)], 


where C is a constant depending on A and the coordinate system. 

When the manifold is sufficiently differentiable every form @ such that 
Aa = 0 is shown to be of class C? and therefore harmonic. The method 
of proof is to represent a in the range of an integral operator whose kernel 


is constructed (see reference 4) with the aid of the parametrix Lr" 7. By 


iteration the kernel can be made square mtegrable 
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For compact manifolds we use a method of Rellich’s’ to prove directly, 
without recourse to a representation as an integral operator, that the 
inverse of A is completely continuous. 

Consider the three linear subspaces of p-forms consisting of, respectively, 
the harmonic forms, the range of d and the range of 6. The self-adjointness 
of dé + éd easily implies that the closures of these three spaces span the 
Hilbert space of p-forms. ‘Together with the other results this leads to the 
standard decomposition theorem (method of orthogonal projections) for 
compact manifolds and a variation thereof for non-compact manifolds. 
It is interesting that a simple theorem of Hilbert space theory (that for a 
self-adjoint operator the closures of the range and the set mapped into zero 
span the Hilbert space) appears as a tool for proving global theorems. 


* Paper done under contract with the Office of Naval Research, N6ori 20, T.O. 25, 
NR 048-083; Mathematics Department, University of Chicago. To be submitted in 
partial fulfillment of the requirements for the Ph.D. degree at the University of Chicago. 

+ The subject here presented was suggested to me in outline by Prof. M. H. Stone. 
In addition, I received valuable aid from Prof. kK. O. Friedrichs. 

' Hodge, W. V. D., “A Dirichlet Problem for Harmonic Functionals, with Applications 
to Analytic Varieties,’’ Proc. London Math. Soc., 36, ser. 2, 72-95 (1932). 

? Kodaira, K., “(Uber die harmonischen Tensorfelder in Riemannschen Mannig- 
faltigkeiten,”’ Proc. Imp. Acad. Tokyo, 20, 186-198, 257-261, 353-358 (1944). 

§ Kodaira, K., ‘“Harmonic Fields in Riemannian Manifolds (Generalized Potential 
Theory),” Ann. Math., 50, 587-665 (1949). 

‘ Bidal, P., and de Rham, G., ‘Les formes différentielles harmoniques,’’ Comment, 
Math. Helv., 19, 1-49 (1946). 

Rham, G., and Kodaira, k., ‘Harmonic Integrals, 
stitute for Advanced Study (1950). 

® We use ‘‘self-adjoint” in the Hilbert space sense, which in the real case means that 
the operator is not only symnietric but that no extension of it is symmetric. 

’ Friedrichs, K., Differential Operators in Hilbert Spaces,” Am. J. Math., 61, 
523-544 (1939) 

* Friedrichs, K., ‘The Identity of Weak and Strong Extensions of Differential Opera- 
tors, Trans. Am. Math. Soc., §5, 182-151 (1944). 

® Rellich, F., Fin Sate uber mittlere Konvergens, Gottingen Nachrichten, 1930, pp. 
30-35. 
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ON THE ANALYSIS OF REPRESENTATIONS OF THE LINEAR 
GROUP 
By Francis D. MURNAGHAN 
INsTITUTO TECNOLOGICO DE AERONAUTICA, RiO DE JANEIRO 
Communicated October 16, 1950 


There is associated with each partition of any positive integer m into not 
more than m parts an irreducible, homogeneous, rational integral repre- 
sentation of degree m of the full linear group of dimension m and all such 
representations of the full linear group are furnished by this association. 
If a is a typical matrix of the linear group of dimension n we denote by 
A = a(A) the matrix which represents a in the irreducible representation 
which is associated with the partition (A) = (Ai, Ae, . ~~, Ax) of m, AL 2 
ho 2... 2 Ne > O,R <n, and we denote by {A} the trace of A. {A} is 
a polynomial function of the power sums 5), So, . . . , Sm of the characteristic 
numbers of a which is known (in honor of I. Schur) as an S-function and 
which is said to be of weight m. If (u) is any partition of any positive 
integer p into not more than .V parts (where N is the dimension of the 
various matrices A) the irreducible representation of the full linear group 
of dimension N which is associated with (u) induces a representation, in 
general reducible, of the group of matrices A and we denote the characters 
of this representation by the symbol {A} @ {uj}. Thus {A} @ {yu} is the S- 
function of the power sums 5S, S:, ..., 5, of the characteristic numbers of 
the matrices A = a(\) which is associated with the partition (yu) of p; it 
is, accordingly, a polynomial function of weight mp of 5, 2, ..., Sm and, 
hence, a linear combination, with non-negative integral coefficients, of S- 
functions of weight mp. The “new product’”’ (so called by Littlewood!) 
{A} ®@ fu} of S-functions is of central importance in the theory of algebraic 
invariants. Littlewood! has given some three methods for carrying out 
this new multiplication of S-functions but he confesses that, even in the 
case where the partition (4) contains only one element, “many difficulties 
arise and it has not so far been possible to overcome them.” In the present 
note we announce some results obtained in studying this problem. 

In the first place we remark that {A} @ {u} may be read off from the 
character table of the symmetric group on mp letters. For example, 


{2} = (S,? + So) where 210) (s,° 55189" + 455), So = 


~ 


20 — + 450). Si? is the square of {31°} which is available in the 


literature’ and .S, is furnished by three columns of the character table of the 
symmetric group on 10 letters (note that we only require a small portion 
of the character table and that this is the portion which is relatively easy to 


2 
f 
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calculate). We obtain S;? = {627} + {6217} + {614} + {532} + [5317] 
4 215221) + 245214) + {519} + f422} + + 264321) + + 
+ + + 13222) + +. {3214}, = (622) — 
16212} + {614} — {582} + {5312} — {515} + {422} — {4212} + [423] — 
— + {37212 {3214} so that {312} @ f2} = 
1627) + 1614} + {5312} + + {5213} + (422) + {4321} + {4313} 
+ {427} + {4214} + [37217]. At the same time we obtain {317} @ 


{12} = — So) = {6212} + {532} + {521} + (5219) + {515} + 


(4212) + {4321} + {4315} + {42°12} + {3°22} 4+ {3214} which is, as it 
should be, the sum of the conjugates of the terms in {317} @ {2}. This 
method is quite practicable for all partitions (A) of m and (4) of p where 
mp < 12. For higher values of mp it becomes impracticable since the 
character table of the symmetric group on r letters is an unwieldy object 
ifr > 12 (it being a 231 X 231 table ifr = 16 and a 627 X 627 table ifr = 
20). 

We denote by {m} the partition of m which contains only one element 
4, = m. Then the formula {m} @ {2} = {2m} + {2m — 2,2} +...+ 
{m?}, if mis even, and {m} @ {2} = {2m} + {2m — 2,2} +...+ {m+ 
1, m — 1}, if m is odd, has been furnished by Littlewood.! ya {m} 
{m} = + {2m — 1,1} +... + [m, m} this implies that {m} @ 
{12} = {2m — 1,1} + (2m — 3,3} +...+ {m+ 1, m — 1}, if miseven, 
and {m} @ {17} = {2m — 1,1} + {2m — 3,3} +... + {m?}, if m is 
odd. ‘To these we add the following two results which are valid for any 
positive integer m > 2: 


I(a). {m — 1, 1} @ {2} = 212m — 2j, 2j} + —j-2,j+ 1,1} 
+ — 2, 2j, 2} 2m — 1, L’}. 

j J 
Each of the four summations starts out with 7 = 1 and continues as long 
as the partitions are in natural (non-ascending) order. Thus {51} @ 


{2} = {10.2} + {921} + {91°} {S4} + {S31} + {827} + {741} + 
(73127) + + {G51} + (642) + {5°17}. 


1(b). {m — 1,1} ® {17} = + 
+ 2j — 3,27 + 1,2} + 
— 2j — 2, 2j, 1%}. 


Thus {61} @ {17} = {10.17} + (93) + (921} + {831} + {82127} + {75} + 
ITAL} + + {651} + (6412) 4 {522}. 


On adding the formulas and we obtain 


| 
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| 
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(2m — 2,14} + 
Lam -j- 2} + 
Thus {61}? = (122) + + + + + 
+ 2410.31] + 110.2") + 195} + 21941) + (982) + (9312) + 
186} + + + {842} + + {74} 24761) + (752) + (751°) + 
+ 
{m — 2,2 2 @ {2) = (2m — 4,24) + — — 
2 2j 4, 2j + 2,2) + — 23-3, 49+ 1,14 + 
2j — 6, 2,4) + — 5,23 +1, 3,1) + 
im — {(m — 1) 2}. 


= {2m — 4, 3, 1} + f2m — 5,3, 2) — 
2-3, +3} + Liem 


2 2j — 5,2 + 3,2 2) + 4 2j + 


2,2) + 2j — 4, 25 + 5,5 + 
2,3} + + 23 — 7,23 + 


1, 2? 


I1(c). im — 2, 2}? = {2m — 4, 3, 1} + {2m — 4,2 2) + 
2 {2m — 9, 3, (2m 5.2%, 1) + 


J 
SU 5,j +: 3,2) + + 
3, 3} + 2 + 


om 3,4) + +23, 1} + 
j 
j 
If the integer m of which (A) is a partition is even and if we denote by 
(A*) the partition of m which is conjugate to, or associated with, (A) then 
{a*} @ | p} is obtained by starring all the terms which appear in the 
analysis of |Aj @ {u} while, if m is odd, {A*} @ {#} is obtained by starring all 


— 
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the terms which appear in the analysis of {A} @ {u*}. Thus the formulas 
we have given serve to determine {A} ® {u} when (u) is (2) or (1?) and (A) 
is any partition of 3, 4 or 5 with the single exception of the self-associated 
partition (31?) of 5. Since we have given above {31°} @ {2} and {31°} @ {17} 
we are in possession of {A} @ {2} and {A} @ {1°} when (A) is any partition 
of 3, 4 or 5 and, generally, if (A) = (m — 1,1), (m — 2,2), (21"~ *), (221" ~ 4). 

The method employed by Littlewood and his coworkers’ * depends on 
the following fact. If {A} is any S-function the derivative of {A} with 
respect to 5; is the sum of the S-functions obtained by subtracting 1, in 
turn, from each element of (A). We denote this derivative of {A} with 
respect to s; by {A}; and it follows at once that [{A} @ fu}]i = [{A} @ 
{u}i] {A}i. From a knowledge of the derivative of {A} @ {u} with respect 
to Littlewood tries to obtain @ but since {A} involves, 
in addition to 5), the variables sy, s3,..., this passage cannot be made with 
certainty. In Litthewood’s words “‘at times a choice of certain alternatives 
is presented, Choosing the correct partitions the work develops rapidly 
to the correct conclusion. But if a wrong partition is chosen the working 
goes further and further astray until the error is apparent . . . much time 
can be lost in following false tracks.'’ These false tracks may be avoided 
when we observe that the relation [}A} @ {uf} = @ {A}. is only 
the first of a series of such relations. If ¢ is any positive integer the deriva- 
tive of {A} with respect to s,/q is the sum of the S-functions obtained by 


subtracting g in turn from each element of (A) (it being understood that 
(A) is extended by the addition of an indefinite number of zeros and that 
disordered partitions are re-ordered in the usual way). We denote this 
derivative of !\} with respect to s,/q by the symbol {A}, and we have 
the following formulas for [{A}@ {u}],,g = 2,3,... 


[1A} @ tuple = @ + @ tafe) @ — 
@ = (LAY @ fads] Ads + @ @ 
+ 
@ = @ Ada + @ — + 


and soon. For example, 


[{2} @ = [{2} @ @ = (2) @ [21] + 
{{2} @ [{2} — — 
= —{2}[{3} — {21} + {17}]; 
{2} ® [{2} — 
O; [{2} = —{2}; 
0,g > 6. 


bo 


to 


bo te 


In conclusion we remark that @ {m} = }5}2M,..., the sum- 
mation being over all partitions (A) of m; {2} @ {m — 1, L{ is obtained by 


| 
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multiplying {2} @ {m — 1} by {2} and subtracting {2} ® {m}. For the 
deeper partitions () of m it seems simplest to use the method outlined above. 
Thus to evaluate {2} @ {m — 2, 2} we suppose that {2} @ {u} has been 
evaluated for all partitions (u) of any integer < m and use the relations 


@ tm — 2, 23h, = {42} [{m — 3, 2} + — 2, {1}, 
[{2} tm — 2, = {2} @ [tm — 3, 2} + {m — 2, 1}] + 
[{m — 4,2} + {m — 2}]} [{2} — {17}] 
and so on.‘ 
' Littlewood, D. E., Phil. Trans. Roy. Soc., Ser. A., 239, 305-365 (1944). 
? Murnaghan, F. D., Am. J. Math., 60, 44-65 (1938). 
3 Foulkes, H.O., J. Lond. Math. Soc., 25, 205-209 (1950). 
“A paper by G. de B. Robinson, treating the problem discussed in the present note, 
has appeared in the current issue of the Canadian Journal of Mathematics. 


THE CHARACTERS OF THE SYMMETRIC GROUP 
By Francis D. MURNAGHAN 
INSTITUTO TECNOLOGICO DE AERONAUTICA, SAO José pos CaAMpos, BRAZIL 


Communicated November 25, 1950 


If (A) = (An, SRL 2 ... > O, is any 
partition of any positive integer m there is associated with (A) an irreducible 
representation of the symmetric group S,,, i.e., the group of m! permuta- 
tions of m letters. The elements of S,, which have the same cycle struc- 
ture, i.e., which contain a unary cycles, a: binary cycles, ..., am cycles on 
m letters, belong to the same class (a) = (a, a, ..., @») and the character 
of any representation of S,, is a class function. The characters of the 
irreducible representations of S,, (each of which is associated with a parti- 
tion of m) are known as the characters of S,, and Frobenius furnished for- 
mulas giving the character of {A} = {Aj, Av, ..., Ax}, ie., of the irreducible 
representation of S,, which is associated with the partition (A) of m, asa 
function of (a) = (a, ..., a») for certain partitions of m. It seems to have 
escaped his attention, and that of subsequent writers on the subject, that 
these and other similar formulas can be procured in a purely mechanical 
tuamner involving little more than the trouble of writing them down; it 
is the object of this note to describe this mechanical process. An account 
of the theory undérlying the process will appear in Lectures on Matrices and 
Matrix Groups to be published by the Centro de Pesquizas Fisicas, Rio de 
Janeiro. 

We first define a process of differentiation of the symbols {A} = {Aj, ... 
x}. If pis any positive integer we denote by {\j, ..., Ax}p the sum 


| 
: 
F 
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where any partition which is disordered by subtraction of p from one of its 
elements is supposed reordered in the usual way; we term {Aj, ..., Ax}p the 
derivative of {A;, ..., Ax} with respect to the integer p. Secondly, if we 
have any polynomial function of a, ..., am We suppose it expressed in 


terms of the binomial coefficients a», and 


we say that its integral with respect to a» is the function obtained by re- 
garding a, a3, ... as constants and replacing any constant C by Cag, any 


constant times a, by that constant times Py any constant times . 


rs and so on. Thus the integral of a, with re- 


by that constant times ( 


spect to a, instead of being 5 02", iS 5 ax(a. — 1) and soon. Similarly the 


integral of any constant with respect to aj is that constant times a3; the 
integral of any constant times a; with respect to a; is that constant times 

=) and so on (and aj, a», a4, ... are treated as constants when we are 
integrating with respect to as). 

After these preliminaries we are ready to explain how to construct the 
formula furnishing the character of {m — p, mi, ..., wei} Where (uw, ..., 
Me-1) = isa partition of p. This formula does not involve apy, -.., 
and the part of it which is independent of ag, ..., a, is obtained by setting 
m = a in Frobenius’ formula for the dimension of {m — p, mw, ..-, Mer}. 
We denote the functions of a, ..., a, which we wish to determine by [y]. 
Thus [0] = 1, since {m} is the identity representation, and [1] = a; — 1 
since the dimension of {m — 1, 1} ism — 1. The part of [2] which is 
independent of ag is youl — 8), since the dimension of {m — 2,2} is vm 
(m — 3) and soon. ‘To obtain the part of [u] which involves a» we differ- 
entiate {u} with respect to 2, replace each {...} in the derivative by the 
corresponding [...] and integrate the result with respect to a: Thus 
{2}. = {0} — [0] = 1 and so the part of [2] which involves az is a. Thus 
the character of {m — 2,2} is furnished by the formula 


[2] = 3) + ae, 


Similarly, since {1°}, = — {0} — —1, the character of {m — 2,1°} is fur- 
nished by the formula 


; 
| 
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For p = 3 we determine the part of [u] which involves a; in the same way; 
we calculate {u}3, replace each {...} in the result by the corresponding 
[...], and integrate the result with respect to a3. If we obtain in this in- 
tegral with respect to a3 a term which already appeared in the integral with 
respect to a2 we do not repeat the term. Thus the part of [3] which is inde- 


pendent of a, and a is peu — 1)(a, — 5), since the dimension of {m — 


3,3} is 5m(m — 1)(m — 5); {3}. = {1} > a — 1 and so the part of [3] 


which involves ay is (a; — l)ae; {3}; = {0} — 1 and so the part of [3] 
which involves a; (and which has not already appeared in the part of [3] 
which involves a2) is a3. Thus the character of the irreducible representa- 
tion {m — 3,3} of S,, is furnished by the formula 


[3] = — 1)(ay — 5) + (ay — + ay. 


The formulas which we derive in this simple mechanical manner have a 
certain universal appeal: They are valid for all symmetric groups Sp, i.e., 
for every value of m. While there is no difficulty in deriving the formulas, 
they take up a great deal of space when p is larger than, say, 6. But this is 
to be expected; the symmetric group is a complex object when m is as 
large as, say, 16. Sig contains 16! elements and has 231 irreducible repre- 
sentations and 231 classes. Thus the mere printing of its character table 
is no light matter, this table being a 231 X 231 table involving (231)? 
entries. We append the formulas furnishing the characters of the irre- 
ducible representations {m — 5, (u)|, where («) is any partition of 5, the 
formulas for p = 3 and p = 4 being already available.' 


1 
[5] = 1) (ay 2)(ay 3) (ay 9) + 1)(ay 5) ay + 


(ay »(%) + + 3) ag + (ay 1) ay + 


| 
[4 1] = — 1)(a, — 2)(ay — 4)(a, — 8) + ~ 2)(ay — — 


+ — 1)(a, — 2)a3 — as; 


[3 2] = — 1)(a, — 2)(a; — 5)(a, — 7) + — 1)(a, — 5)ay + 


(a, — »(%) + aa, — — 1)(a, — 2)a3 — (ay — 1)ay; 


1 
[3 37} 1)(ay ngs 3) (ay 4) (ay = 7) + (ay = 1) ay 
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— (2) 


1 
(2? 1] = 1) (ay 3) (ay 5) (ay 6) 1)(ay 2) (ay 


3) ae + (ay »(%) 3) as + (ay 1) a4; 


[2 14] = 2)(ay 3) (a4 4) (ay 6) 2)(a, 4) ay 


+ — — as; 


[1°] = 129°" 1)(a 2) (a4 3) (ay 4) (ay 5) 1) (ay 


2)(a1 3) a + (ay »(%) — 1)(ay 2) a3 


(ay > 1) + 5. 
1 Murnaghan, F. D., Theory of Group Representations, The Johns Hopkins Press (1938). 


VECTOR FIELDS ON THE n-SPHERE 
By N. E. Steenrop J. H. C. Wuirenrap 
DEPARTMENT OF MATHEMATICS, PRINCETON UNIVERSITY AND OXFORD UNIVERSITY 
Communicated by S. Lefschetz, November 18, 1950 


1. Introduction.—Let S" denote the n-sphere, and let 2* be the highest 

power of 2 dividing m + 1: 
n+ 1 = 2*(2r + 1). (1) 
Our main result is 

THEOREM 1.1. Any set of 2* continuous vector fields tangent to S” are 
somewhere dependent. 

Note first that if r = 0 in (1) then 1.1 is trivially true. 

The case k = 0 of 1.1 is the classical result that a vector field on a sphere 
of even dimension has at least one zero, The case k = 1 was proved by 
Eckmann! and also by G. Whitehead.*?_ A proof was given for the case k = 
2 by G. Whitehead,* but was based on the erroneous assertion of Pon- 
trjagin* that = 0. 

Whether S" admits 2* — 1 independent vector fields is unsettled. It is 
known": ? that, fork = 1, 2 and 3, S" admits 2* — 1 such fields if m + 1 = 
0 mod 2". These fields are constructed using the complex, quaternionic 
and Cayley division algebras, respectively. It is not known whether there 
exist real division algebras on 2* units for any k = 4. 
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Before proving 1.1, we state various consequences. As shown by one of 
us,® 1.1 is equivalent to 

THEOREM 1.2. Jf n and k areas in (1), and if 2* ¢ q Sn — 2, then S" 
does not admit a continuous field of tangent q-dimensional hyperplanes. Nor 
does S" admit a continuous quadratic form (symmetric covariant tensor of 
order 2) which, at each point, is non-singular and of signature q (1.e., q nega- 
tive characteristic roots). 

Let C"*' be cartesian (n + 1)-space of coordinates (4, ..., tn+1); let S" 
be the locus }77"' ¢,2 = 1; let C" be the subspace of the first m variables of 
(i.e., tars =... = bar, = 0); let = C™ NA S*; let R,,; denote 
the rotation group (orthogonal transformations of determinant +1) of 
C"*t!; let R, be the rotation group of C” regarded as the subgroup of Ryy1 
which leaves pointwise fixed the (n — m + 1)-space orthogonal to C”; and 
let x, denote the vector (or point), all of whose components are zero save 
t, = 1. Define the projection 


p: — by a@e (2) 


Then, for x €.S", p~'(x) is a left coset of R,; so p is just the fibering of R,,1 
by left cosets of R,, and S” is the left coset space R,,;/R,. We define the 
Stiefel manifold V,,,, ,,, to be the left coset space 


V Rave (q < n). 3) 


Using inclusion relations among cosets, decomposes into the composition 
of two natural projections 


g h 


An ordered, orthonormal set of g + 1 vectors in C"*! will be called a (q + 
1)-frame. R,,, is transitive on the set of all (¢ + 1)-frames based at the 
origin, and R,,-, leaves fixed the frame (%, g.1, Xavi). Thus Vass, gut 
can be interpreted as the manifold of (g + 1)-frames based at the origin of 
C"*!, The projection / assigns to each (¢ + 1)-frame the end-point of its 
last vector. ‘To each (gq + 1)-frame corresponds the g-frame tangent to S" 
obtained by a parallel translation of the first g vectors to the end-point of 
the last vector. Using this correspondence, we may interpret Vays, a1 as 
the manifold of all q-frames tangent to S", and h as the assignment of the 
point of tangency. By the usual process of orthogonalization, q fields of 
tangent vectors, which are independent at each point, can be reduced to q 
fields which are orthonormal at each point, i.e., to a field of tangent q- 
frames. Thus 1.1 is equivalent to 

THEOREM 1.3. If n and k are as in (1) with r > O, then the fiber bundle 
h: Vasa, 28s: — S” does not have a cross-section (i.e., there does not exist a 
continuous map f: S" such that hf(x) = x for all x). 
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If f: S*—> R,4 is a cross-section for R,,; — S", then gf is a cross-section 
for Vass. hence 

Coroiiary 1.4. If n is not of the form 2* — 1, then the fiber bundle 
— does not have a cross-section, 

Using the result of Kirchhoff,’ we have 

Corouvary 1.5. Jf n is not of the form 2* — 2, then S" does not admit a 
complex analytic structure. Even more, it does not admit an exterior form of 
degree 2 which is non-singular at each point. 

As shown by one of us,* 1.4 implies 

Coro.iary 1.6, [fis not of the form 2* — 1 and B is an n-sphere bundle, 
then the image of m,(S") in m,(B) under the natural map is cyclic of even or 
infinite order. Thus the fiber S" is not contractible to a point in B. 

Coroiary 1.7. If nis not of the form 2* — 1,andm > n, then S” cannot 
be fibered by n-spheres. 

1.8. A representation of S” as an n-sphere bundle over 
is possible only if n has the form 2* — Landm = 2**! — J, 

In the last three corollaries, the group of the bundle is R,,,. It is not 
known if they are valid with R,,, replaced by the group of all homeomor- 
phisms of S”. 

2. Squaring Operations in Projective Spaces. Before proving 1.1, we 
shall establish a preliminary result. Let P" be real projective n-space. 
It is known that the cohomology group //“(P") with coefficients mod 2 is a 
cyclic group of order 2 with generator u(0 S q Sn). The cohomology 
ring 11(P") = /1*(P") satisfies = For any complex K, let 
Sq’: — H**(K) be the squaring homomorphism determined by 
cup-(q — j) products.* 

2.1. Jn P", we have Sq/u" = c(j, q)u"*! whereO <j < 


< 
q+ j Sn, and c(j, q) ts the mod 2 value of the binomial coefficient e! 


The case q = | follows from the general properties :’ = and 
=ununifg = dinu. Using the general formula of H. Cartan? 


Sq/(u uv) = (Sq'u) u 
and the proposition’ Sq‘u = Ofori > dim u, we proceed by induction: 
= Sq/(u' 1)" = + Sq'u' yu 


To apply 2.1, we shall develop a simple criterion for evaluating c(j, q). 
It follows from 2.1 that c(j, g) is the mod 2 coefficient of w in the expansion 
of (1 + u)* where 1 = uw°andu = u'. Let 


i 
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be the dyadic expansions of j, q (a,b, = Oor 1). Calculating mod 2, we 
have (1 + u)? = 1 + w?; hence (1 + u)* = 1 + u*"; and therefore 


(1 + u)* = + (5) 
By the uniqueness of the dyadic expansion, the coefficient of u/ in (5) is 
bi... bi, where 4, ..., 1, are those values of i for which a, = 1. Therefore 
LeMMA 2.2. c(j, q) = 1 af and only tf the non-zero terms of the dyadic ex- 


pansion of j form a subcollection of those of q. 
THEOREM 2.3. Ifn +1 = 2*(2r + 1) with r > O, then, in P", we have 


= u", and = 0 for 0 <j < 2. 


To prove this, note first that 
mn — = — = 4 — 1) = 2! + — J), 


Therefore 2* is a non-zero term of the dyadic expansion of nm — 2*. Then 
2.2 implies that ¢(2*, nm — 2*) = 1, and the first assertion follows from 2.1, 
To prove the second, observe that 0 < 7 < 2* implies that j has the form 
2°(2t-+ 1) whereO Ss <k. Since 


n= — 1 + = 2! + 


it follows that the coefficient of 2° in the dyadic expansion of n — j is 0, 
while it is 1 in that of 7. Then 2.2 implies c(j, m — 7) = 0; and 2.1 applies 
to complete the proof. 

3. Proof of 1.1.—-We adopt the notations preceding 1.3, and assume that 
r>0O. Let £”, E™ denote the hemispheres of S” defined by t,,,., 2 0, 
tmvi & O, respectively. Their totality, form = 0, 1, ..., ", forms a cell 
complex on S". 

For each x € S", let T(x) € R,,, be the composition of reflection in the n- 
plane through the origin perpendicular to x, followed by reflection in the n- 
plane through the origin perpendicular to x. Clearly 7(x) = 7T(—x). 
Since x,,; is fixed under the first reflection, p7(x) is the image of x,,; 
under the second reflection. It follows that p7 maps each great circle 
through x,,; onto itself with a covering map of degree 2... From this we 
obtain the following properties of 7: 


T(x) = T(x’) if and only if x = x’ orx = ~—x’, (6) 
T(S") = P" isa real projective n-space in R,,,,, (7) 
T(S") = = P*n R,,,; form = 0,1, ..., (8) 


pT(S"") = and pT maps — (9) 
topologically onto — 


‘ 
= | 
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Let T’ = T\S""'. By (9), we may regard 7” as a map S"~! — R,. 
The maps 7, 7’ were first introduced by Hopf," and 7” is called the char- 
acteristic or canonical map for the bundle — G. Whitehead'’ 
has shown that S” admits g continuous, independent vector fields if and 
only if 7’ is homotopic in R, toa map of S"~'into R,,. Let 


be the natural map, i.e., g’ = gR,. Since g’ is a bundle projection, 7” is 
deformable into a map if and only if g’T’ is homotopic to a 
constant. Taking g = 2‘, we have only to show that g’7’: — V,, 
is not homotopic to a constant. 

One of us!* has given a representation of V,, 2* as a cell complex A con- 
sisting of one 0-cell e® = no m-cells for 0 < m <n — 2%, exactly 
one m-celle” = g/T’(E”) forn — 2* Sm Sn — 1, and exactly one cell for 
each set of distinct integers m, ..., m, between n — 2* — 1 and n, whose 
dimension is }oim,. The cell of least dimension greater than m —1 corre- 
sponds tom, = n and m, = n — 2*+ 1. It has dimension 


2n — + 1 = n + 2*(2r — 1), see (1). 


Asr > 0, it follows that the n-skeleton AK” of A coincides with K"~'. As 
g’T’ is a cellular map, we have g’T7’(S"~!) = K"~'. Since A” = K"~'!, we 
have that g’7” is homotopic to a constant in V’,, 9 if and only if it is homo- 
topic to a constant in A"~'. Thus we have only to prove that g’7”, re- 
garded as a map S”'-» AK"~'!, is not homotopic to a constant. 

Let e” = gT(E?) © Vass. 1, and let K’ = e" uy Let gi: P"— K’ 
be the map g P". If we regard P” as a cell complex composed of the cells 
T(E"),m = 0,1, ...,, then g, is cellular and for n — S m S n maps the 
open m-cell 7(/%) topologically onto the open m-celle”. It follows that g, 
induces isomorphisms of the mod 2 cohomology groups 


git: H'(K’) HP") for n—-2 sisn. (10) 
Since g,*Sq’ = Sq’gi*, the result of 2.3 applies with (10) to give 
Sq w"-* = w" (11) 


where w’' is the non-zero element of J/'(4K’’). 

Suppose, as is impossible, that g’7’: S"~!—> A"~! is homotopic to a con- 
stant. This means that g7 regarded as a map /. — K’ can be extended to 
a map f: S" K’ such that f(/") ¢ Clearly f*w” is the non-zero 
element of 1/"(S"); and ftw"-“ = 0 since 7"-“(S") = 0. But Sq/f* = 
f*Sq’; so we have a contradiction with (11), and the proof is complete. 

ReMARK |. If, in the preceding proof, we replace 2* by 7 < 2*, then (11) 
becomes Sq/w"? = 0 (see 2.3), and we obtain no contradiction. This 
shows that the bound 2‘ in 1.1 is the least value provided by our argument. 
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REMARK 2. Theorem 1.1, with S" replaced by real projective space P", 
has been proved by Stiefel!‘ and by Hopf. Their result is a corollary of 
1.1; for the double covering S" —> P” enables one to “‘lift’’ vectors on P” 
into vectors on .S” preserving independence. It may seem surprising that 
the theorem for P® should be easier than for S". However the proof for P" 
uses the fact that a certain characteristic class is non-zero. When 7 is odd, 
all the characteristic classes of S" are zero. Thus the primary obstructions 
vanish, and we have shown that an obstruction of higher order is non-zero. 
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MANY-FIGURE APPROXIMATIONS TO 2, AND 
DISTRIBUTION OF DIGITS IN AND 


By Horace S. UHLER 
YALe UNIVERSITY 
Communicated by J. B. Whitehead, November 15, 1950 


The present investigation was begun with the single intention of testing 
a value for 4/2 as computed by J. Marcus Boorman about 65 years ago. 
This datum was chosen because the author had acquired the impression 
that Boorman seidom if ever produced a reliable result in his calculations 
and it seemed fair to check this particular approximation since it was 
claimed that he had verified many of his digits. More precisely we read! 
“Mr. J. M. Boorman of New York has computed (pp. 112-115) the square 
root of 2 to 520 decimals, and verified 486 of them.’ Then the 486 decimal 
digits are quoted. 

At the time the author was unaware of the recent work of René Coustal? 
in which the approximation was carried to 1032 guaranteed decimal places. 


- i 
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This noteworthy contribution stimulated the checking of Coustal’s answer 
and the extension of it to at least 1542 decimal places by the present critic. 

The work was begun by squaring Boorman’s result to find out to what 
extent if any the square differed from 2. This operation was conveniently 
effected by treating his approximation as the sum of two numbers, 6 and 
c. 6 represented the first 243 digits of B’s 1/2 with the decimal point 
temporarily disregarded. ¢ comprised the last 243 figures of which the 
first figure (8) occupied the 244th decimal place. (The 243rd_ place 
happens to be filled by zero.) The-right-hand digit of b? occupied the 
484th decimal place while the respective leading figures of 2bc and c? fell 
in the 243rd and 487th places. This calculation showed that (b + c)? 
consisted of 1, the decimal point, and 315 nines followed by the weird 
sequence 

59999 99987 99999 99999 99999 99999 99999 99960 00000 

08000 08000 00000 00000 00000 00000 00000 00000 00000 

00000 00069 99999 99998 80000 00000 00000 59999 99999 

99999 99999 99999 99996 43244 95623 60074 24492 20770 

92230 99222 52571 61930 68368 12237 83658 55185 33491 

64379 =19306 83618 72765 59353 71418 55794 54954 
Hence it has now been proved that Boorman’s revised value of \/2 was 
correct to 315 decimals but erroneous for the remaining 171 figures. 

The next stage of the present work comprised the following steps, the 
theory and justification for which will be outlined below. The preceding 
value of (b + c)? was subtracted from 2, this difference was divided by 4, 
and then [2 — (6 + c)?]/4 was multiplied by the 315-place correct approxi- 
mation to \/2. In this way the author obtained the value of \/2 to 611 
decimal places. 

About the same time Dr. John W. Wrench, Jr., computed an approxi- 
mation to 1/2. On August 27, 1950, he juxtaposed his value against the 
611-place work strip and found complete agreement as far as his last digit 
in the 460th decimal place. 

Let a denote a definite approximation to / N and ¢ the correction to a 
to give a higher approximation. More precisely 


VN =a+te (1) 


(Expressed in decimals \/N and ¢ would never terminate when N is not a 
perfect square.) In practice we square a and subtract the result from N. 
Let this difference be symbolized by 6 so that 


N—a’® = 6. (2) 
Now « = VN — a from equation (1) and YN = a/Va*/N = 
a(l — 6/N)~“* from equation (2). The condition for convergence 
6°/N*? < 1 should be amply fulfilled so that the required practical expansion 
may be written 


| 
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ab { 35 35 
€= on an ant (3) 


In the present case N = 2 and the first order term applied earlier is ¢ = 
a(6/4). 


The author's 612 digit approximation was identical with the correspond- 


TABLE 1 
V2 = 
1.41421 35623 73095 04880 16887 24209 69807 85696-71875 87604 
80731 76679 73799 07324 78462 10703 88503 87534 32764 15727 
35013 84623 09122 97024 92483 60558 50737 21264 41214 97099 
93583 14132 22665 92750 55927 55799 95050 11527 82060 57147 
01095 59971 60597 02745 34596 86201 47285 17418 64088 91986 
09552 «32923 «0484308714. 3214508397 52514-07989 
68725 33965 46331 '80882 96406 20615 25835 23950 54745-75028 
77599 61729 838557 52203) 87531 85701 13543 74603 40849 88471 
60386 89997 06990 04815 03054 40277 90316 45424 78230 68492 
we— 
93691 86215 80578 46311 15966 68713 01301 56185 68987 23728 
52885 09264 86124 94977 15421 83342 04285 68606 01468 24720 
77143 58548 74155 65706 96776 53720 22648 54470 15858 = 80162 
07584 74922 65722 60020 855 46652 14583 98893 94437 09265 
91800 31138 82464 68157 08263 01005 94858 70400 31864 80342 
19489 72782 90641 04507 26368 81318 73985 52561 17322 04024 
50912. 27700» 22694-11275) 73627) 2804957381) (08967) 50401 83698 
68368 45072 57993 64729 06076 29969 41380 47565 48237 28997 
18032 68024 74420 62926 91248 59052 18100 44598 42150 59112 
02494 41341 72853 14781 05803 60337 10773 09182 86931 47101 
71111 68391 65817 26889 41975 87165 82152 12822 95184 88472 


Ce 

08969 46338 «62891 56288 27659 52685 | 22676 =53239 69461 
75112. 91602. 40871-55101) 85150) 45588) 12875) 60052-63146 80171 
27402 65396 «94702 40300-51749) 53188) 63138) 51881 63478 
00156 93691 76881 85237 86840 52287 83762 93892 14300 65586 
95686 85964 59515 55016 50983-68960 36887 82311 438804 
15576 65104 08839 14292 33811 32060 52433 62948 53170 49915 
77175 62285 49741 43899 91880 21762 43096 52065 64211 82731 
67262 57539 72559 34637) 23863 22614 2742 62220 86711 
55839 59992 65211 76252 69891 75409 88159 34864 00884 57085 
18147 22318 14204 07042 65090 56532 33339-84364 «57865 79679 
65192 67292 39987 53666 17215 98257 88602 63363 6178 


ing portion of Coustal’s longer value. This agreement justified the 
application of formulae (2) and (3) to the 612 figure value of a. The 
second calculation was carried to the 1040th decimal place. , The result 
previously obtained by Coustal agreed absolutely with this slightly more 
extended datum to 1034 places. This circumstance fully confirms Coustal’s 
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claim that ‘Tous ces chiffres sont exacts, la 1033° décimale est sirement 
un zéro, la 1034° est probablement un 5." Hence ‘‘probablement’’ may 
now be changed to “certainement.” 

In this case 6 = (0.161859---) &K = (0.572260---) X 
1.e., 611 decimally antecedent zeros; 3a6?/32 = (3/32)(1.414213---) 
(1.618---)* = (0.347346---) X 10-12, 1222 ante-zeros 
which safely include the 1034th decimal place. 

Finally in order to contribute something entirely new, while using a 
machine in the same non-electronic class as that presumably employed 
by Coustal, the perfect 1033 digits of \/2 were assigned toa. Then the 
1544 decimal approximation to ./2 as given below was derived from 


TABLE 2 

Py Pn 
240 2068 
1804 0544 
2125 0059 
3831 0470 
0853 0039 
1634 0176 
3539 0527 
3703 0091 
4945 0468 
4833 0394 
5196 0846 
4980 . 2254 
5504 0993 
5638 0407 
5869 0056 


woe we wan 


= 


x 


x 


equations (2) and (3). The superior letters B, W, and C in table 1 mark 
the limits of accuracy attained, respectively, by Boorman, Wrench and 
Coustal. 

In addition to checking throughout all stages of the work with small auxili- 
ary moduli of the form 10" + 1 the final answer was carefully squared. This 
long multiplication gave an unbroken succession of 1544 nines between 
the decimal point (after 1) and a 2. The writer guarantees 1542 decimal 
figures and he believes that the digits in the 1543rd and 1544th places are 
correct. 

In the last case 6 = (0.153376---) & 10~'8*; a6/4 = (0.542267---) X 
10-5, ie., 1033 decimally antecedent zeros; = (3/32) 
(1.414218: - -)(1.533764- ++)? & 107% = (0.311891---) 107%, 
2066 ante-zeros which overwhelm the 1544th decimal place. 


” X_* 

100 

200 
300 
400 

500 

600 

700 

800 

900 

1000 8 

1100 
1200 a 316 
1300 
1400 
1500 9.12 

1545 9 
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Table 2, given above, relative to the distribution of the digits of \/2 and 
1/\/2 may possess statistical interest. The values of p, were obtained 
by linear interpolation in a table of percentage points of the x? distribution 
by Catherine M. Thompson.* However the value of P for nm = 1 and nine 
degrees of freedom could not be located by the writer, hence to avoid 
extrapolation he calculated the product 


P= 

with the aid of a thirteen-place table of logarithms by H. Andoyer (1922), 
and found p = 0.000562497302. 

In table 2 the second and third columns pertain to \/2 and they indicate 
settling down to a fairly steady normal state for n > 800. The fourth 
and fifth columns refer to 1/ 1/2 and, with two or three exceptions, they 
show noteworthy small values of x,? and p,. The relatively large values of 
pn for n = 1100, 1200 and 1300 can be traced back to very appreciable 
steady increases in the number of 4’s combined with simultaneous de- 
ficiencies in the number of 9’s occurring. 

1 Editor, Mathematical Magazine, 1, 164 (1887). 


2 Coustal, René, Compt. rend., 230, 431, 432 (1950). 
3 Thompson, Catherine M., Biometrika, 32, 187-191 (1941-1942). 
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